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ABSTRACT 


The  first  thermal  imagery  from  the  Spatially  Enhanced  Broadband  Array 
Spectrograph  System  (SEBASS)  was  analyzed  for  target  detection  purposes.  Data 
were  acquired  during  Exercise  WESTERN  RAINBOW  in  October  of  1995  with 
Aerospace  Corporation’s  new  instrument.  SEBASS  is  a  thermal  infrared  spectral 
imaging  system  which  collects  data  in  128  wavelength  bands  from  7.8  to  13.4  pm. 
The  system  has  128  spatial  pixels  and  a  ground  sample  distance  of  2.5  feet  under 
the  operating  conditions  of  the  exercise.  The  data  reduction  process  included 
thermal  calibration  using  on-board  calibration  sources,  atmospheric  correction 
through  the  use  of  in-scene  data,  and  derivation  of  thermal  emissivity.  Analysis  of 
the  data  was  performed  using  techniques  adopted  from  hyperspectral  imagery.  The 
application  of  the  principal  components  transformation,  the  spectral  angle  mapper, 
and  a  spectral  matched  filter  showed  substantial  enhancement  of  target  detection 
capabilities  when  compared  to  the  analysis  of  standard  infrared  imagery. 
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1.  INTRODUCTION 


Rapid  advances  in  remote  sensing  technologies  have  led  to  new  sensors  which 
greatly  improve  our  ability  to  gather  information  about  distant  objects  without  direct 
physical  contact.  While  many  of  these  sensors  and  technologies  were  designed  for 
environmental  monitoring  and  earth  science,  they  may  have  important  uses  in  the 
Department  of  Defense  (DOD)  as  today’s  warfighter  strives  to  gain  a  complete  picture  of 
the  battlespace  around  him.  As  the  numbers  and  capabilities  of  these  sensors  continue  to 
grow,  the  evaluation  of  their  capabilities  in  a  support  to  military  operations  (SMO)  setting 
has  become  increasingly  important. 

Hyperspectral  imagery  is  one  of  the  most  promising  of  these  emerging 
technologies.  Hyperspectral  imaging  refers  to  a  version  of  imaging  spectroscopy  in 
which  measurements  of  reflected  or  emitted  light  are  taken  in  hundreds  of  narrow  spectral 
bands  (Multispectral  Users  Guide,  1995).  Because  the  reflection  or  emission  of 
electromagnetic  energy  is  characteristic  of  a  material’s  molecular  makeup,  this  high 
resolution  trace  of  radiation  versus  wavelength  can  form  a  graphical  record  unique  to  a 
given  material  (Rinker,  1990).  The  high  spectral  resolution  of  hyperspectral  sensors  can 
provide  opportunities  to  detect  subtle  differences  between  material  signatures.  From  the 
standpoint  of  military  targeting  and  identification,  hyperspectral  sensors  have  a  potential 
beyond  many  previous  remote  sensing  instruments. 

Most  hyperspectral  instruments,  like  the  Airborne  Visible/Infrared  Imaging 
Spectrometer  (AVIRIS)  and  the  Hyperspectral  Digital  Imagery  Collection  Experiment 
(HYDICE),  collect  information  in  the  region  of  the  electromagnetic  spectrum  dominated 
by  solar  illumination.  These  instruments  acquire  a  complete  reflectance  spectrum  by 
recording  the  reflected  solar  radiation  from  an  object.  Although  these  sensors  have 
achieved  success  in  material  identification,  they  are  ineffective  at  night  in  the  absence  of 
solar  illumination.  The  Aerospace  Corporation  created  the  Spatially  Enhanced 
Broadband  Array  Spectrograph  System  (SEBASS)  to  extend  the  range  of  hyperspectral 
sensors  into  the  thermal  infrared  region.  In  contrast  to  the  visible  region  of  the 
electromagnetic  spectrum,  the  infrared  region  is  dominated  by  thermal  self-emission. 
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This  allows  the  sensor  to  operate  equally  well  during  both  day  and  night.  Thermal 
infrared  hyperspectral  sensors  have  the  potential  to  perform  target  detection,  terrain 
categorization,  and  materials  classification  regardless  of  illumination  conditions.  They 
achieve  this  by  examining  the  emissivity  differences  between  target  objects  in  the  same 
way  that  visible  hyperspectral  sensors  examine  differences  in  solar  reflectivity.  Because 
of  its  freedom  from  the  restraint  of  solar  illumination,  SEBASS  may  represent  a  step 
forward  in  the  SMO  capabilities  of  hyperspectral  sensors. 

The  analysis  of  hyperspectral  imagery  differs  dramatically  from  the  literal 
interpretation  of  single-band  imagery  because  of  the  nature  of  hyperspectral  data.  By 
recording  measurements  in  hundreds  of  spectral  bands,  hyperspectral  sensors  generate 
large  amounts  of  data  that  require  complex  methods  of  analysis  (Multispectral  Users 
Guide,  1995).  Unlike  single-band  systems  in  which  target  detection  is  based  on  contrast 
and  spatial  resolution,  hyperspectral  sensors  discriminate  between  various  target  and 
background  types  on  the  basis  of  their  reflective  or  emissive  characteristics.  Many 
techniques  have  been  devised  to  accomplish  this  spectral  analysis.  The  overall  shape  of  a 
target’s  radiation  curve  can  be  compared  to  known  laboratory  measurements.  The 
locations  of  specific  features  within  several  spectral  curves  may  be  recorded  and 
compared.  Even  statistical  analysis  of  the  variations  between  all  of  the  spectral  data 
within  a  entire  scene  can  lead  to  significant  information  about  scene  content.  The  choice 
and  application  of  a  particular  technique  has  a  direct  impact  on  the  ability  of  the 
sensor/algorithm  combinaition  to  detect  or  classify  a  target. 

The  MASINT  Hyperspectral  Study  (MHS)  was  initiated  in  1993  to  assess  the 
utility  of  hyperspectral  imagery  (MHS:  Interim  Progress  Report,  1995).  The  MHS 
defined  three  SMO  goals  to  guide  its  efforts: 

•  Find  threat  vehicles  and  defeat  Camouflage,  Concealment,  and  Deception 
(CCD)  efforts  better  than  currently  done 

•  Improve  current  capabilities  for  the  categorization  of  terrain  features 

•  Demonstrate  the  utility  of  thermal  hyperspectral  products  for  military  and 
terrain  categorization  problems  (MHS:  Interim  Progress  Report,  1995). 
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In  support  of  these  goals,  the  MHS  program  sponsored  Exercise  WESTERN 
RAINBOW  at  the  Yuma  Proving  Ground,  Arizona,  in  October  1995.  During  this 
exercise,  a  variety  of  airborne  and  ground-based  hyperspectral  sensors  were  employed  in 
a  SMO  scenario. 

This  thesis  will  present  an  assessment  of  SEBASS  and  its  capability  to  conduct 
SMO  missions.  It  evaluates  the  utility  of  thermal  hyperspectral  products  for  military  and 
terrain  categorization  problems  using  data  collected  during  WESTERN  RAINBOW.  A 
description  of  hyperspectral  remote  sensing  theory  and  hyperspectral  data  analysis 
techniques  is  presented  in  Chapter  H.  Chapter  m  is  a  discussion  of  the  history  and 
development  of  SEBASS.  The  experimental  setup  of  Exercise  WESTERN  RAINBOW 
is  outlined  in  Chapter  IV  while  the  experimental  observations  collected  during  the 
exercise  are  described  in  Chapter  V.  Chapter  V  also  includes  ground  truth  measurements, 
SEBASS  data  set  descriptions,  and  a  complete  description  of  the  SEBASS  data 
preparation  process.  Spatial  and  spectral  analysis  of  the  WESTERN  RAINBOW  data  is 
discussed  in  Chapter  VI.  Finally,  Chapter  VTI  will  present  the  conclusions  drawn  from 
the  analysis  of  SEBASS  data  collected  during  the  exercise. 
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II.  INFRARED  HYPERSPECTRAL  REMOTE  SENSING 


The  goal  of  spectral  remote  sensing  is  the  identification  of  material  through  the 
examination  of  reflected  or  emitted  electromagnetic  radiation.  All  materials  have  a 
unique  spectral  signature  based  on  their  ability  to  reflect  or  emit  electromagnetic  energy 
at  different  wavelengths.  The  spectral  signature  of  an  object  is  usually  presented  as  a 
reflectance  or  emittance  spectrum.  A  reflectance  spectrum  graphically  depicts  the  ability 
of  an  object  to  reflect  incident  electromagnetic  radiation  at  different  wavelengths.  An 
emittance  spectrum  depicts  the  ability  of  a  material  to  emit  thermal  energy. 

Figure  2.1  depicts  the  differences  between  the  reflectance  spectrum  of  healthy 
vegetation  and  two  types  of  camouflage  (Multispectral  Users  Guide,  1995).  The  figure 
illustrates  how  spectral  information  can  be  used  to  identify  and  discriminate  between 
materials  based  on  their  spectral  properties.  From  0.4  to  0.7  pm,  all  three  materials 
appear  similar.  This  wavelength  region  corresponds  to  the  visible  portion  of  the 
electromagnetic  spectrum  and  confirms  that  modem  camouflage  is  designed  to  simulate 
natural  vegetation.  At  wavelengths  longer  than  0.7  pm,  significant  differences  between 
the  man-made  camouflage  and  natural  vegetation  begin  to  become  discernible.  Sensors 
able  to  record  electromagnetic  radiation  at  these  longer  wavelengths  would  be  able  to 
discriminate  between  all  three  types  of  objects. 

The  ability  of  a  sensor  to  detect  subtle  spectral  differences  between  objects 
depends  on  several  factors.  As  shown  in  the  discussion  of  Figure  2.1,  different  regions 
of  the  electromagnetic  spectmm  can  reveal  differences  between  materials  that  may  not 
otherwise  be  detectable.  Also,  a  high  spectral  resolution  sensor  will  be  able  to  detect 
small,  unique  features  within  a  particular  spectmm.  For  example,  the  sensor  used  to 
collect  the  data  for  Figure  2.1  has  a  spectral  resolution  better  than  0.05  pm.  This  narrow 
bandwidth  allows  the  sensor  to  accurately  measure  the  narrow  drop  in  reflectivity  at  1.45 
pm  in  the  spectmm  of  green  vegetation.  The  atmosphere  between  an  object  and  the 
sensor  can  also  affect  the  quality  of  the  reflectance  or  emittance  spectmm  collected  by  an 
instmment.  Molecules  within  the  atmosphere  absorb,  emit,  and  scatter  electromagnetic 
energy.  Sensors  which  are  designed  to  exploit  the  differences  between  spectra  collected 
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from  ground  materials  must  have  a  method  of  compensating  for  these  degrading  effects 
which  can  mask  the  true  emission  characteristics  of  the  target  material.  Finally,  spectral 
data  analysis  techniques  vary  in  their  ability  to  detect  and  exploit  features  within  collected 
spectral  data. 

This  chapter  presents  a  description  of  infrared  hyperspectral  remote  sensing.  It 
begins  with  a  description  of  the  principles  of  spectral  remote  sensing.  The  specific 
characteristics  of  infrared  remote  sensing  are  then  presented.  Finally,  it  concludes  with  a 
discussion  of  several  classes  of  hyperspectral  data  analysis  techniques. 

A.  SPECTRAL  REMOTE  SENSING 

Spectral  sensors  are  designed  to  exploit  particular  regions  of  the  electromagnetic 
spectrum  in  an  effort  to  take  advantage  of  the  properties  of  electromagnetic  radiation  at 
different  wavelengths.  A  sensor  may  be  designed  to  operate  during  periods  of  low  solar 
illumination,  to  detect  a  characteristic  spectral  feature  of  a  target  material,  or  to 
discriminate  between  several  types  of  minerals.  To  satisfy  these  different  uses,  the 
number  of  spectral  bands  measured  by  remote  sensing  instruments  as  well  as  their  width 
and  location  within  the  electromagnetic  spectrum  will  vary.  These  characteristics  are 
selected  to  achieve  the  level  of  resolution  required  to  satisfy  the  sensor’s  intended  goal. 

1.  Electromagnetic  Spectrum 

The  portion  of  the  electromagnetic  spectrum  used  most  frequently  in  spectral 
remote  sensing  ranges  from  0.2  pm  to  14.0  pm.  Figure  2.2  graphically  depicts  this  region 
of  the  spectrum  (Multispectral  Users  Guide,  1995).  Sensors  operating  in  the  region  from 
0.2  pm  to  2.5  pm  collect  and  record  reflected  solar  energy.  In  the  transition  region 
between  2.5  pm  and  8.0  pm,  both  reflected  solar  energy  and  thermal  self-emission 
contribute  to  the  electromagnetic  energy  measured  by  a  sensor.  Between  8.0  and  14.0 
pm,  thermal  self-emission  dominates  the  observed  spectrum.  Sensors  operating  at 
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wavelengths  longer  than  8.0  fim  collect  emitted  radiation  and  do  not  depend  on  reflected 
solar  illumination. 

a.  Visible,  Near  Infrared,  and  Shortwave  Infrared 

Spectral  imagery  in  the  visible,  near  infrared  (NBR),  and  shortwave 
infrared  (SWIR)  region  of  the  electromagnetic  spectmm  is  dominated  by  the  reflection  of 
solar  energy.  Visible  light  ranges  from  0.4  ^im  to  0.7  ^m  and  is  usually  classified 
according  to  color.  The  shorter  wavelengths  correspond  to  blue  light  while  the  longer 
wavelengths  correspond  to  red  light.  The  reflectance  of  many  materials  differ  within  the 
visible  region  providing  the  many  variations  of  color  we  observe  between  objects 
(Anderson,  et  al,  1994).  The  NIR  (0.7  to  1.1  p,m)  and  SWIR  (1.1  to  2.5  ^im)  regions  of 
the  electromagnetic  spectrum,  although  not  detectable  by  the  human  eye,  are  also 
dominated  by  the  reflection  of  solar  energy. 

Sensors  operating  in  the  visible,  NIR,  and  SWIR  regions  exploit  variations 
in  the  ability  of  objects  to  reflect  radiation.  Energy  in  the  blue  region  of  visible  light  is 
absorbed  by  chlorophyll.  It  also  can  penetrate  clear  water  to  a  depth  of  about  40m 
(Multispectral  Users  Guide,  1995).  Reflected  energy  in  the  visible  red  region  has  been 
useful  for  vegetation  and  soil  discrimination  (Multispectral  Users  Guide,  1995).  NIR 
reflectance  is  used  for  vegetation  analysis  because  it  is  strongly  affected  by  the  cellular 
structure  of  leaf  tissue.  Sensors  operating  in  the  SWIR  region  are  used  to  detect  moisture 
content  and  highlight  the  contrast  between  vegetation  types. 

b.  Midwave  Infrared  and  Longwave  Infrared 

Thermal  emission  of  radiation  becomes  a  significant  factor  in  the  midwave 
infrared  (MWIR)  and  longwave  infrared  (LWIR)  regions  of  the  electromagnetic 
spectrum.  MWIR  radiation  includes  wavelengths  from  3.0  to  5.0  p,m.  The  LWIR  region 
extends  from  8.0  to  14.0  |xm.  In  these  longer  wavelength  regions,  materials  emit  more 
radiation  than  they  reflect  from  the  sun.  As  a  result,  sensors  which  operate  within  this 
region  are  not  constrained  by  solar  illumination.  Variations  in  the  ability  of  materials  to 
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emit  radiation  are  the  thermal  infrared  counterpart  to  reflectance  variations  that  are  used 
for  material  discrimination  in  visible,  NIR,  and  SWIR  remote  sensing. 

During  the  day,  MWIR  energy  is  a  mix  of  reflected  and  emitted  energy. 
At  night,  the  reflected  solar  component  is  reduced  and  the  energy  recorded  by  a  sensor 
operating  in  the  MWIR  region  consists  almost  entirely  of  emitted  radiation.  MWIR 
imagery  can  be  used  to  penetrate  smoke  and  has  been  used  to  measure  surface 
temperature  (Multispectral  Users  Guide,  1995).  LWIR  imagery  has  traditionally  been 
used  for  thermal  temperature  measurements.  LWIR  sensors  and  emittance  spectra  can 
also  be  used  to  conduct  vegetation  classification  and  perform  geologic  mapping 
(Multispectral  Users  Guide,  1995). 

2.  Spectral  Resolution 

The  spectral  resolution  of  a  sensor  describes  its  ability  to  detect  and  record 
differences  in  reflected  or  emitted  electromagnetic  energy  as  a  function  of  wavelength. 
Spectral  resolution  varies  from  sensor  to  sensor  based  on  intended  use  and  current 
technology.  Low  resolution  sensors  include  single-band  sensors  which  collect  and  record 
all  electromagnetic  energy  within  a  broad  spectral  band  as  well  as  sensors  which  collect 
data  in  several  wide  wavelength  bands.  High  resolution  sensors  record  hundreds  and 
even  thousands  of  energy  measurements  in  very  narrow  bands. 

Single-band  remote  sensing  instruments  record  the  total  electromagnetic  energy 
received  within  a  single  wavelength  band.  Panchromatic  photography  and  forward- 
looking  infrared  (FLIR)  sensors  are  examples  of  single-band  sensors  (Anderson,  et  al, 
1994).  Panchromatic  photography  is  sensitive  to  visible  wavelength  energy  between  0.4 
p.m  and  0.7  )Lim.  FLIR  sensors  generally  operate  between  8|iim  and  12  |im.  Bright  and 
dark  areas  within  these  images  indicate  high  and  low  surface  reflectances  or  warm  and 
cold  surfaces. 

The  drawback  of  single-band  sensors  is  that  many  materials  will  appear  similar 
within  a  single,  wide  band.  The  spectral  features  that  distinguish  one  material  from 
another  may  be  very  narrow  or  may  be  located  in  widely  dispersed  regions  of  the 
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electromagnetic  spectrum.  To  increase  the  possibility  of  detecting  these  spectral 
differences,  sensors  which  record  the  same  scene  in  separate  wavelength  bands  are  used. 
Multi-band  sensors  are  classified  according  to  the  number  and  width  of  their  spectral 
bands.  The  three  classes  of  multi-band  sensors  are  multispectral,  hyperspectral,  and 
ultraspectral.  Figure  2.3  depicts  these  three  classes  of  sensors  and  their  relative 
capabilities  (Multispectral  Users  Guide,  1995). 

a.  Multispectral  Sensors 

Multispectral  sensors  have  2-20  spectral  bands  (Anderson,  et  al,  1994). 
These  sensors  take  multiple  images  of  a  scene  by  collecting  and  recording  radiation 
within  each  separate  band.  The  bands  selected  for  multispectral  sensors  do  not  need  to  be 
contiguous  and  often  include  a  combination  of  visible,  NIR,  SWIR,  MWIR  and  LWIR 
spectral  bands.  By  comparing  the  recorded  electromagnetic  energy  within  these  bands, 
multispectral  sensors  can  perform  terrain  classification,  camouflage  detection,  soil 
analysis,  and  trafficability  analysis  (Anderson,  et  al,  1994).  A  variety  of  airborne 
multispectral  sensors  have  been  developed  including  the  current  generation  of  Landsat 
and  SPOT  sensors.  Table  2.1  is  a  summary  of  the  bands  employed  in  the  Landsat 
Thematic  Mapper  multispectral  instrument. 


1  Band  Number  Spectral  Bands  (mm)  I 

1 

0.45  -  0.52  (blue) 

2 

0.52  -  0.60  (green) 

3 

0.63  -  0.69  (red) 

4 

0.76  -  0.90  (NIR) 

5 

1.55  - 1.75  (SWIR) 

6 

10.4-12.5  (LWIR) 

7 

2.08  -  2.35  (SWIR) 

Table  2.1.  Landsat  Thematic  Mapper  Spectral  Bands.  From  Richards,  1993. 


b.  Hyperspectral  Sensors 

Hyperspectral  sensors  record  data  in  hundreds  of  narrow,  contiguous 
spectral  bands  (Multispectral  Users  Guide,  1995).  "The  increased  number  of  sensor 
bands  provides  higher  spectral  resolution  and  more  opportunities  to  detect  subtle  spectral 
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differences  in  signatures  that  are  too  narrow  to  be  differentiated  on  multispectral 
imagery.”  (Multispectral  Users  Guide,  1995)  This  increased  spectral  resolution  allows 
the  instrument  to  approximate  the  true  reflectance  or  emittance  spectrum  of  target  objects 
and  allow  for  the  identification  and  classification  of  surface  materials  on  a  pixel-by-pixel 
basis  (Goetz,  et  al,  1985).  NASA's  AVIRIS  sensor  (Vane,  et  al  ,1993),  the  DOD 
HYDICE  instrument  (Rickard,  et  al,  1993),  and  Aerospace's  SEBASS  (Hackwell,  et  al, 
1996)  are  examples  of  airborne  hyperspectral  sensors.  AVIRIS  and  HYDICE  operate  in 
the  visible,  NIR,  and  SWIR  regions  of  the  electromagnetic  spectrum.  SEBASS  collects 
data  in  the  LWIR  region  of  the  spectrum.  The  basic  characteristics  of  these  sensors  are 
listed  in  Table  2.2. 


Instrument 

Spectral  Range  (um) 

Number  of  Spectral  Bands 

AVIRIS 

0.4  -  2.5 

224 

HYDICE 

0.4  -  2.5 

206 

SEBASS 

7.8-13.4 

128 

Table  2.2.  Spectral  Band  Characteristics  of  Several  Hyperspectral  Sensors.  After  Vane, 
et  al,  1993;  Rickard,  et  al,  1993;  and  Hackwell,  et  al,  1996. 

c.  Ultraspectral  Sensors 

Ultraspectral  sensors  record  data  in  thousand  of  contiguous  spectral  bands. 
They  are  useful  for  detecting  narrow  spectral  phenomena  within  a  small  portion  of  the 
electromagnetic  spectmm.  This  very  high  resolution  capability  has  been  used  for  specific 
material  identification  and  molecular  identification  of  aerosol  and  plume  components 
(Multispectral  Users  Guide,  1995). 

B.  INFRARED  REMOTE  SENSING 

A  spectral  sensor  collects  reflected  and  emitted  radiation  from  a  source.  The 
initial  energy  source  may  be  the  sun,  thermal  emission  from  the  target  object,  or  energy 
radiated  from  the  sensor.  In  the  LWIR  region  of  the  spectrum,  the  source  of  radiation  is 
the  target  material  itself.  The  radiation  recorded  by  the  sensor  is  the  result  of  interactions 
between  the  source  energy  and  the  atmospheric  constituents  between  the  sensor  and  the 
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target.  The  goal  of  infrared  remote  sensing  is  to  determine  characteristics  of  the  source 
object  by  working  backwards  from  the  measured  radiance  at  the  sensor.  This  requires  an 
understanding  of  the  emissive  characteristics  of  materials  as  well  as  the  radiative  transfer 
process  by  which  energy  propagates  through  the  atmosphere. 

1.  The  Planck  Function 

All  materials  with  a  temperature  above  absolute  zero  emit  radiation.  To  describe 
the  manner  in  which  an  object  radiates  energy,  the  concept  of  a  perfect  radiation  emitter, 
or  blackbody,  has  been  adopted.  The  amount  of  radiation  emitted  by  a  blackbody  is  a 
function  of  the  temperature  of  the  blackbody  and  the  wavelength  of  emitted  energy. 
Figure  2.4  is  a  graph  of  the  radiation  characteristics  of  several  blackbodies  at  various 
temperatures.  As  the  temperature  of  a  blackbody  increases,  the  amount  of  enndtted 
radiation  also  increases.  The  figure  also  shows  that  the  radiation  emitted  from  a 
blackbody  is  not  constant  across  the  electromagnetic  spectrum,  but  rises  to  a  maximum 
and  then  gradually  declines.  The  wavelength  at  which  this  peak  of  emitted  radiation 
occurs  decreases  as  the  temperature  of  the  blackbody  increases.  The  equation  governing 
blackbody  radiation  is  called  the  Planck  Radiation  Law  and  is  presented  as  Equation  2.1. 


Equation  2.1  jr] 

e^-l 

C,  =  1.19104X  10'®  \y.W ! cm^\im  ^r] 

C2  =  1.438769x10^  k] 

X,=  wavelength  [pn] 

T = temperature  [^] 


The  emission  characteristics  of  real  materials  are  described  by  comparing  the 
thermal  emission  of  a  real  object  to  a  blackbody  radiating  at  the  same  temperature.  The 
ratio  of  these  quantities  is  the  spectral  emissivity  of  an  object  as  defined  in  Equation  2.2. 
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Equation  2.2  e,  = 

’■  B^(X) 

4«to<rd  =  emitted  radiance 

The  spectral  emissivity  of  a  material  also  describes  its  ability  to  absorb 
electromagnetic  energy  according  to  Kirchoff’s  law  which  states,  "..for  equilibrium 
conditions  the  absorptivity  of  a  sample  is  exactly  equal  to  its  emissivity.  "(Zissis,  1993) 

2.  Molecular  Properties 

Electronic,  vibrational,  and  rotational  transitions  within  atoms  and  molecules 
account  for  the  absorption,  emission,  and  reflection  properties  of  materials.  Spectral 
features  in  the  visible  and  ultraviolet  regions  of  the  electromagnetic  spectrum  are  the 
products  of  electronic  processes.  Within  the  infrared  region,  spectral  features  are  the 
results  of  molecular  vibrations  and  rotations. 

In  general,  atomic  gases  produce  discrete  absorption  and  emission  features.  Table 
2.3  shows  the  prominent  absorption  features  of  several  important  gases.  The  discrete 
nature  of  absorption  by  gases  produces  the  distinctive  features  in  atmospheric 
transmission  described  below. 


Species 

Wavelengthliim) 

CO 

4.66 

CO2 

2.7,  4.2, 14.99 

N2O 

4.5,7.78,16.98 

H2O 

2.73,6.27,15 

O3 

9.01,9.59,14.18 

Table  2.3.  Prominent  Absorption  Features  of  Atmospheric  Gases. 
After  Anderson,  et  al,  1994,  and  Zissis,  1993. 


Liquids  and  solids  produce  broad  absorption  features.  The  absorption  features  of 
solids  and  liquids  within  the  0.4  -  14.0  pm  region  of  the  electromagnetic  spectrum  are 
due  to  electronic  and  vibrational  transitions  at  resonant  frequencies.  Between  0.4  and  2.5 
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fxm,  spectral  features  are  generally  caused  by  electronic  energy  transitions  within  atoms 
and  molecules.  Above  2.5  ^im,  vibrational  absorption  resonances  are  the  primary  source 
of  spectral  features.  For  a  more  complete  description  of  the  optical  properties  of 
materials,  see  Anderson,  et  al  (1994). 

Liquid  water  has  a  high  absorption  coefficient  and  is  essentially  opaque  above  2.5 
M'tn.  Figure  2.5  shows  the  reflectance  of  water  for  the  2.0  to  14.0  |j.m  region  (Anderson, 
et  al,  1994).  Within  this  region,  the  reflectance  of  water  is  very  low,  ranging  from  0.5% 
to  4.0%.  This  low  reflectance  indicates  that  the  majority  of  incident  radiation  is  absorbed. 
A  water  layer  as  thin  as  0. 1  mm  is  enough  to  absorb  almost  all  longwave  radiation.  As  a 
result,  even  a  thin  layer  of  moisture  on  a  target  can  significantly  affect  the  amount  of 
radiation  emitted  or  reflected  from  it. 

Silicate  minerals  exhibit  strong  spectral  features  in  the  LWIR  region.  The  nature 
of  Si-O  bonds  causes  fundamental  vibrational  absorption  resonances  above  8.0  |iim. 

These  resonances,  called  reststrahlen  bands,  appear  as  minima  in  emissivity  spectra. 
Figure  2.6  displays  examples  of  emissive  spectra  from  four  different  silicate  minerals 
(Zissis,  1993).  The  reststrahlen  feature  can  be  seen  in  all  four  spectra  as  a  minimum  in 
the  8.0  to  1 1.0  (im  region.  The  location  of  the  reststrahlen  feature  is  different  for  each 
mineral  because  of  varying  percentages  of  silicate.  This  information  can  be  used  to 
classify  and  discriminate  between  various  minerals. 

Vegetation  displays  near-blackbody  qualities  in  the  thermal  infrared  region.  Early 
measurements  indicated  that  vegetation  displayed  very  little,  spectral  information  in  this 
region.  More  accurate  and  extensive  data  has  recently  indicated  that  although  the  spectral 
variations  may  be  small,  significant  spectral  information  is  present  (Anderson,  et  al, 

1994).  The  emissivity  of  vegetation  can  be  altered  by  water  content  and  degree  of 
senescence.  Figure  2.7  is  a  plot  the  reflectance  spectra  of  healthy  and  senescent  grass 
(John  Hopkins  University,  1995).  The  difference  in  overall  shape  between  these  two 
spectra  is  evident. 
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3. 


Radiative  Transfer 


Spectral  sensors  record  the  amount  of  energy  received  from  a  particular  direction. 
They  commonly  measure  incident  energy  in  terms  of  radiance.  Radiance  is  the  amount  of 
power  per  unit  area  that  is  incident  on  a  surface  from  a  particular  direction.  In  the 
infrared  region  of  the  spectrum,  radiance  measured  at  a  sensor  is  a  combination  of  energy 
emitted  by  the  target  object,  energy  emitted  by  the  surrounding  background,  and  energy 
emitted  or  absorbed  by  the  atmosphere  between  the  object  and  the  sensor  (Multispectral 
Users  Guide,  1995).  Figure  2.8  depicts  the  various  paths  along  which  energy  from  these 
sources  energy  may  travel  to  a  sensor.  Path  1  represents  the  target  source  energy.  As 
described  by  Equations  2.1  and  2.2,  the  energy  emitted  by  an  object  is  a  function  of  its 
temperature  and  emissivity.  Not  all  of  the  energy  emitted  by  an  object  reaches  the  sensor. 
A  fraction  is  absorbed  by  atmospheric  constituents  and  some  is  scattered  in  directions 
away  from  the  sensor.  To  account  for  these  effects,  an  additional  factor,  atmospheric 
transmittance,  is  introduced.  Atmospheric  transmittance  is  a  measure  of  how  much 
radiation  will  propagate  through  the  atmosphere  to  the  sensor.  The  contribution  of  energy 
from  path  1  is  described  by  Equation  2.3. 

Equation  2.3  [}XlT/cm^|Xm  sr] 

'C;j^  =  atmospheric  transmittance 

The  second  source  of  energy  reaching  the  sensor  represents  thermal  energy 
emitted  by  the  atmosphere  onto  an  object  and  reflected  toward  the  sensor.  The  energy 
emitted  by  the  atmosphere  toward  the  target  surface  is  called  downwelling  radiance.  The 
percentage  of  downwelling  radiance  that  reaches  the  sensor  is  found  by  multiplying  the 
downwelling  radiance  by  the  reflectivity  of  the  material  and  the  transmittance  of  the 
atmosphere.  Within  the  thermal  region,  the  reflectivity  of  most  materials  is  generally 
small  and  can  be  calculated  by  subtracting  the  emissivity  of  an  object  from  one.  The 
contribution  of  path  2  is  described  by  Equation  2.4. 
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Equation  2.4  L2 ^x)^downweiimg  [|iW/cm^|iun  ^r] 

^dcrnnweiung  =  atmosphedc  downwelling  radiance 

Finally,  energy  emitted  by  the  atmosphere  may  travel  directly  toward  the  sensor  as 
depicted  by  path  3.  This  upwelling  radiance  is  a  function  of  the  atmospheric  constituents 
and  ambient  air  temperature.  The  contribution  of  upwelling  radiance  is  described  in 
Equation  2.5. 

Equation  2.5  L,  =  [\iWlcm^\m  sr] 

^upwelling  =  atmospheric  upwelling  radiance 

Equations  2.3  through  2.5  can  be  combined  into  a  single  equation  representing  all 
of  the  electromagnetic  energy  measured  by  a  sensor.  This  equation  is  the  radiative 
transfer  equation  for  LWIR  remote  sensing  and  is  presented  as  Equation  2.6.  A  more 
complete  discussion  of  infrared  radiative  transfer  is  given  in  Kidder  and  Haar  (1995). 

Equation2.6 

jr] 


4.  Atmospheric  Effects 

As  mentioned  above,  the  atmosphere  affects  the  radiance  measured  at  the  sensor 
in  several  ways.  Atmospheric  constituents  can  attenuate  and  scatter  solar  radiation  in  the 
visible  portion  of  the  spectrum.  In  the  thermal  infrared  region,  the  molecules  in  the 
atmosphere  emit  radiation  and  serve  as  a  diffuse  source  of  energy.  General  path  losses 
from  the  target  to  the  sensor  can  also  alter  the  radiance  measurements  at  the  sensor. 

Figure  2.9  shows  the  transmittance  of  the  atmosphere  throughout  most  of  the 
electromagnetic  spectrum  (Kidder  and  Haar,  1995).  Several  ‘windows’  in  which  a  large 
percentage  of  source  radiation  will  propagate  are  evident.  These  are  largely  defined  by 
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absorption  due  to  water  vapor  and  CO2.  The  8.0  to  14.0  inm  window  is  punctuated  by 
ozone  absorption  at  9.6  pm  and  by  water  vapor  absorption  throughout  the  band. 

The  goal  of  atmospheric  compensation  is  to  accurately  determine  atmospheric 
transmittance  as  well  as  upwelling  and  downwelling  radiance.  Several  methods  have 
been  devised  to  determine  these  parameters.  Atmospheric  compensation  methods  range 
from  the  use  of  detailed  atmospheric  models  to  the  use  of  experimentally  collected  data. 

a.  Atmospheric  Models 

The  cumulative  effects  of  atmospheric  attenuation  and  scattering  can  be 
determined  using  radiative  transfer  models.  Most  atmospheric  models  accept  a  standard 
set  of  input  parameters  that  describe  the  environmental  state  and  produce  output 
parameters  such  as  upwelling  radiance  and  transmittance.  Multispectral  and 
hyperspectral  sensor  applications  conunonly  employ  a  low  resolution  atmospheric 
transmittance  model  (LOWTRAN)  and  a  moderate  resolution  model  (MODTRAN) 
(Anderson,  et  al,  1994;  Zissis,  1993).  To  increase  the  flexibility  of  atmospheric 
modeling,  Aerospace  Corporation  developed  the  Atmospheric  Radiance  and 
Transmittance  Spectra  model  (ARTSPEC)  which  estimates  atmospheric  parameters  in  the 
2.8  - 14.0  ^im  range  (Burke,  1995).  ARTSPEC  gains  its  flexibility  through  the  use  of  a 
large  variety  of  environmental  inputs.  The  environmental  parameters  used  as  inputs  for 
these  models  are  usually  collected  on-site  with  meteorological  instruments  or  estimated 
using  the  acquired  spectral  data. 

b.  Empirical  Line  Method 

The  empirical  line  method  has  been  used  extensively  in  the  visible,  NIR, 
and  SWIR  regions  of  the  spectrum  to  derive  surface  reflectance  (Farrand,  et  al,  1994). 

The  empirical  line  method  is  based  on  the  following  equation: 
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Equation  2.7 


^sensor  =  PxA  +  B 


[|xW/cm^|im  5r] 


=  surface  material  reflectance 
A  =  multiplicative  term  (x,,  and  instrument  gain) 

B  =  additive  term  (path  radiance  and  instrument  offset) 

The  empirical  line  method  solves  for  the  multiplicative  and  additive  terms 
in  each  wavelength  band  by  examining  one  light  and  one  dark  target  within  the  scene. 
The  reflectance  spectra  of  the  two  target  materials  must  be  known.  Target  reflectances 
are  typically  determined  through  laboratory  analysis  or  the  use  of  field  spectrometers. 
Figure  2.10  demonstrates  the  empirical  line  method  for  a  single  band.  The  slope  of  the 
line  drawn  between  the  two  points  is  the  multiplicative  term,  A .  The  vertical  offset 
represents  the  additive  term,  B .  The  process  is  repeated  for  each  band  of  a  multi-band 
imager.  The  resulting  additive  and  multiplicative  terms  for  each  band  are  then  applied  to 
the  entire  data  set  to  derive  reflectance  spectra  for  every  pixel  in  the  scene. 

c.  Plastic  Ruler  Atmospheric  Compensation  Technique 

The  plastic  ruler  atmospheric  compensation  method  was  developed  by 
Hackwell  and  Hayhurst  (1995)  to  derive  atmospheric  parameters  for  infrared 
hyperspectral  collections  using  in-scene  data.  The  method  is  based  on  a  simplified 
version  of  the  infrared  radiative  transfer  equation  (Equation  2.6).  When  the  emissivity  of 
a  material  equals  1.0,  Equation  2.6  simplifies  to  Equation  2.8. 

Equation  2.8  ='^,5,(7')  +  ^upwelling  [p,W/cm^)Xm  jr] 

This  equation  is  similar  to  Equation  2.7  which  forms  the  foundation  of  the 
empirical  line  method.  The  plastic  ruler  atmospheric  compensation  algorithm  attempts  to 
determine  the  slope  ,Xj^ ,  and  offset ,  ,  of  a  set  of  points  representing  in-scene  data. 

The  process  is  complicated  by  the  fact  that  the  amount  of  radiation  emitted  from  each 
data  point  in  the  scene,  (T) ,  is  a  function  temperature.  The  temperature  of  each  point 
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must  be  estimated  or  measured  prior  to  the  solution  of  the  above  equation.  In  addition, 
determining  the  slope  and  offset  of  a  cluster  of  points  extracted  directly  from  collected 
data  is  more  difficult  than  simply  drawing  a  straight  line  because  the  points  do  not  all 
represent  scene  elements  with  emissivities  close  to  1.0. 

Figure  2.1 1  can  be  used  to  explain  the  plastic  ruler  algorithm.  The  figure 
is  a  scatter  plot  of  data  values  in  a  single  wavelength  band  taken  from  a  SEB  ASS 
collection  mission  during  WESTERN  RAINBOW.  For  each  data  point,  the  surface 
temperature  has  been  estimated  and  the  amount  of  radiation  emitted  by  a  blackbody  at 
this  estimated  temperature,  {T) ,  has  been  calculated.  Each  point  is  then  plotted  on  a 

scatter  plot  of  measured  radiance  versus  blackbody  radiance.  The  plastic  ruler  method 
solves  for  the  atmospheric  transmittance  and  upwelling  radiance  by  drawing  a  line 
through  the  upper  bound  of  the  scatter  plot.  The  points  within  the  upper  bound  represent 
scene  elements  with  emissivities  close  to  1.0.  These  points  satisfy  the  simplified 
radiative  transfer  equation  presented  as  Equation  2.8.  The  plastic  ruler  method  employs 
Kolmogorov-Smimov  statistics  and  noise  data  from  the  instrument  to  find  this  upper 
bound. 

The  initial  assumption  that  the  emissivity  of  each  scene  element  equals  1.0 
somewhere  within  its  spectmm  results  in  relative  atmospheric  transmittance  and 
upwelling  radiance  values  instead  of  absolute  corrections  (Hackwell  and  Hayhurst,  1995). 
Burke  (1995)  combined  the  plastic  ruler  method  with  their  ARTSPEC  atmospheric  model 
to  derive  absolute  atmospheric  corrections.  The  results  from  the  plastic  ruler  method 
were  used  to  guide  the  ARTSPEC  atmospheric  model.  The  ARTSPEC  model  input 
parameters  were  varied  until  the  model’s  results  exhibited  the  same  relative  shape  as  the 
plastic  ruler  results.  The  actual  ARTSPEC  model  results  could  then  be  used  for 
atmospheric  compensation  or  as  scaling  factors  for  the  plastic  ruler  results. 

C.  HYPERSPECTRAL  DATA  ANALYSIS 

A  variety  of  data  analysis  techniques  have  been  developed  to  exploit 
hyperspectral  data.  Hyperspectral  analysis  algorithms  have  been  developed  to  find 
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matches  for  a  particular  target  spectra  within  a  scene,  to  detect  anomalous  spectra  within 
a  scene,  and  to  classify  the  materials  within  a  scene  based  on  spectral  properties.  Three 
representative  algorithms  and  techniques  were  applied  to  the  data  collected  by  SEB  ASS 
at  WESTERN  RAINBOW .  These  techniques  included  the  principal  components 
transformation,  the  spectral  angle  mapper,  and  a  spectral  matched  filter. 

1.  Principal  Components  Analysis 

Hyperspectral  image  data  consist  of  a  series  of  images,  each  representing  the 
reflection  or  emission  of  radiation  within  a  particular  frequency  band.  Because  of  the 
narrow,  contiguous  nature  of  hyperspectral  bands,  a  high  degree  of  correlation  exists 
between  these  spectral  images  (Santisteban  and  Munoz,  1977).  The  principal 
components  transformation  is  a  rotational  transformation  that  projects  the  hyperspectral 
data  into  a  new  coordinate  system  which  is  based  on  the  variance  of  the  original  data. 
When  the  transformation  is  applied  to  hyperspectral  data,  the  spectral  information  of  the 
original  image  is  transformed  into  a  new  set  of  orthogonal  bands.  The  result  is  a  series  of 
decorrelated  images,  called  principal  component  images.  The  number  of  images  equals 
the  number  of  original  spectral  bands. 

A  pictorial  representation  of  the  principal  components  transformation  is  presented 
in  Figures  2.12  and  2.13.  This  example  is  taken  from  Richards  (1993).  These  figures 
demonstrate  the  process  on  a  two-band  data  set.  The  two  figures  are  scatter  plots  of 
several  data  points  plotted  as  band  intensity  values.  The  data  in  Figure  2.12  exhibit  a 
high  degree  of  correlation  between  the  bands.  An  increase  in  band  1  is  almost  sufficient 
to  predict  an  increase  in  band  2.  The  goal  of  principal  components  analysis  (PCA)  is  to 
rotate  the  data  into  a  coordinate  system  in  which  the  correlation  between  the  bands  is 
zero.  Figure  2.13  is  an  illustration  of  such  a  rotation. 

The  PCA  technique  is  accomplished  through  the  Karhunen-Loeve  transformation 
(Santisteban  and  Munoz,  1977).  The  transformation  is  a  linear  transformation  which 
generates  a  new  set  of  coordinates  by  weighting  the  original  spectral  bands.  The  first 
principal  component  maximizes  the  variance  of  the  original  data  into  a  new  band  which  is 
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a  linear  combination  of  measurements  in  the  original  spectral  bands.  The  second  and 
subsequent  principal  components  each  maximize  the  remaining  variance  in  the  data.  The 
transformation  is  accomplished  through  Equation  2.9.  The  linear  operator,  G ,  is  defined 
as  the  transposed  matrix  of  eigenvalues  of  the  scene  covariance  matrix.  The  decorrelation 
of  the  data  is  described  by  Equation  2.10.  The  new  covariance  matrix  is  a  diagonal 
matrix  of  the  eigenvalues  of  the  original  covariance  matrix. 

Equation  2.9  JP.  =  GSi 

Pi  =  vector  of  principal  component  values  of  i'*  pixel 
G  =  linear  transformation 
5,  =  spectrum  of  i'*  pixel 

Equation  2.10  C'  =  GCG' 

C  =  scene  covariance  matrix 
G  =  transposed  matrix  of  eigenvectors  of  C 
C'  =  diagonal  matrix  of  eigenvalues  of  C 

The  first  principal  component  represents  the  largest  data  variance  within  the 
scene.  The  successive  components  contain  less  and  less  of  the  scene  variance.  As  a 
result,  the  first  few  principal  components  will  contain  most  of  the  information  in  the 
original  data.  The  last  few  components  will  contain  mostly  noise.  This  type  of 
transformation  has  been  considered  a  data  compression  technique  because  it  concentrates 
the  original  information  contained  in  hundreds  of  bands  of  spectral  data  into  10  or  20 
principal  component  bands. 

Once  the  linear  transformation  matrix,  G ,  has  been  generated,  it  can  be  applied 
quickly  and  efficiently  to  additional  data  sets  through  simple  matrix  multiplication.  This 
efficiency  makes  PCA  attractive  from  an  operational  perspective.  Unfortunately,  the 
transformation  matrix  is  sensitive  to  changes  in  scene  content.  Another  drawback  of  PCA 
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is  that  the  quantification  of  what  each  principal  component  physically  represents  is 
difficult  since  the  transformation  is  based  strictly  on  variability  of  the  data  set. 

2.  Spectral  Angle  Mapper 

The  spectral  angle  mapper  (SAM)  is  a  spectrum  matching  algorithm  that  labels 
scene  elements  based  on  their  similarity  to  a  particular  target  spectrum  (Kruse,  et  al, 
1993).  The  SAM  algorithm  chosen  for  this  analysis  is  implemented  in  the  Environment 
for  the  Visualization  of  Images  (ENVI)  data  analysis  software  package.  The  algorithm 
compares  the  spectral  data  from  an  image  to  selected  target  or  ‘end  member’  spectra  and 
produces  a  set  of  images  in  which  intensity  represents  the  similarity  to  each  target 
spectrum.  The  algorithm  also  produces  a  classification  image  showing  the  best  match  for 
each  pixel  in  a  scene. 

The  SAM  determines  the  similarity  between  spectra  by  calculating  the  “spectral 
angle”  between  them.  It  was  developed  for  spectral  data  in  the  visible,  NER,  and  SWIR 
regions  of  the  electromagnetic  spectrum.  The  spectra  are  treated  as  vectors  in  a  space 
with  dimensionality  equal  to  the  number  of  spectral  bands.  Figure  2. 14  illustrates  the 
concept  of  “spectral  angle”  for  a  two-dimensional  case.  Each  measurement  is  represented 
as  a  vector  from  the  origin  to  a  point  plotted  in  band  intensity  space.  A  vector 
terminating  close  to  the  origin  represents  a  point  with  little  illumination  while  a  vector 
terminating  far  from  the  origin  represents  a  brightly  illuminated  point.  In  either  case,  the 
direction  of  the  vector  is  characteristic  of  the  material’s  spectral  properties. 

The  SAM  uses  the  angular  separation  between  a  target  spectrum  and  a  reference 
spectral  vector  to  determine  spectral  similarity  (Schowengerdt,  1994).  This  allows  the 
SAM  to  perform  equally  well  regardless  of  illumination  and  gain  conditions.  The  angular 
separation  between  an  unknown  spectrum  and  the  reference  spectrum  can  be  calculated 
by  Equation  2.11. 
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Equation  2.11  a,-  =  cos 

a,-  =  spectral  angle  of  z'*  pixel 
T  =  reference  spectrum 
Sj  =  spectrum  of  V’'  pixel 

The  output  of  the  SAM  is  a  new  image  set  with  the  number  of  images  equal  to  the 
number  of  target  reference  spectra.  Each  image  consists  of  the  spectral  angle  separation 
between  the  pixels  in  the  scene  and  a  particular  target  spectral  vector.  Bright  areas 
indicate  close  matches.  The  use  of  a  minimum  separation  angle  can  be  employed  to 
classify  spectra  as  a  match  to  a  particular  reference  spectra. 

The  SAM  algorithm  can  be  applied  directly  to  radiance  measurements  in  the 
reflective  region  of  the  electromagnetic  spectrum.  In  the  thermal  infrared  region,  the  data 
must  be  converted  to  emissivity  or  apparent  emissivity  prior  to  the  application  of  the 
SAM  algorithm.  This  is  because  the  effects  of  temperature  variations  are  not  simply 
additive  as  in  the  case  of  illumination  variations  in  visible,  NIR,  and  SWIR  data. 

3.  Spectral  Matched  Filter 

The  spectral  matched  filter  (SMF)  is  another  spectral  matching  technique.  It  was 
adapted  from  matched  filter  techniques  used  in  the  spatial  and  temporal  domains.  The 
SMF  technique  used  for  this  analysis  was  implemented  in  the  ENVI  data  analysis 
software  package.  As  described  by  Billingsley,  et  al  (1995),  the  SMF  algorithm  is  a  two- 
stage  process.  The  first  step  characterizes  the  dominant  spectral  components  in  the  scene 
and  ‘whitens’  them.  The  second  step  is  a  search  for  a  match  to  a  spectral  reference  vector 
in  the  residual  subspace. 

The  SMF  uses  scene  spectral  statistics  and  a  spectral  reference  vector  to  calculate 
weighting  factors  which  combine  the  appropriate  spectral  bands  to  attain  the  highest 
signal-to-clutter  ratio  (Stocker,  et  al,  1990).  The  SMF  algorithm  is  implemented  through 
Equation  2.12. 
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Equation  2.12  yi  =b‘  C"’ x, 

y,  =  filtered  intensity  of  i'*  pixel 

b  =  spectral  reference  vector 
C  =  background  spectral  covariance  matrix 
Xi  =  spectram  of  i"'  pixel 


The  SMF  can  also  be  viewed  as  a  weighted  combination  of  principal  components 
as  described  by  Stocker,  et  al  (1990).  This  interpretation  is  presented  in  Figures  2.15  - 
2. 17.  Figure  2.15  illustrates  a  Gaussian  clutter  density  and  an  additive  spectral  signal,  b  . 
Correlation  between  the  two  spectral  bands  causes  the  clutter  density  to  be  elongated 
along  a  major  axis.  Projections  of  the  data  beside  each  waveband  axis  demonstrate  that 
the  signal-to-clutter  ratio  in  either  band  is  not  sufficient  to  detect  the  signal.  A  principal 
component  transformation  results  in  Figure  2.16.  The  transformation  has  enabled  the 
spectral  signal,  b  ,  to  become  distinct  from  the  clutter  in  the  second  principal 
component.  The  SMF  gains  a  further  increase  in  the  signal-to-clutter  ratio  by  forming  the 
optimum  weighted  combination  of  principal  components  using  a  weight  vector  based  on 
the  reference  spectrum.  The  result  is  illustrated  in  Figure  2.17. 

The  output  of  the  SMF  is  a  two-dimensional  spatial  image.  Each  point  in  the 
image  is  the  result  of  an  optimum  weighted  combination  of  spectral  band  measurements. 
High  intensity  pixels  within  the  filtered  image  indicate  the  presence  of  the  reference 
spectra.  A  thresholding  algorithm  can  be  used  to  identify  likely  targets. 
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III.  THE  DEVELOPMENT  OF  THE  SPATIALLY  ENHANCED  BROADBAND 
ARRAY  SPECTROGRAPH  SYSTEM  (SEBASS) 

Thermal  remote  sensing  development  has  followed  a  progression  toward  higher 
spectral  resolution  in  the  drive  toward  identifying  materials  based  on  their  emissive 
spectra.  The  thermal  infrared  multispectral  scanner  (TIMS)  demonstrated  that 
contiguous,  moderately  narrow  bands  could  be  used  to  discriminate  between  various 
minerals  based  on  their  spectral  features  (Richards,  1993).  Bongiovi,  et  al  (1995) 
demonstrated  that  hyperspectral  LWBR  measurements  could  detect  fine  spectral  features 
and  identify  specific  materials  based  on  their  spectral  characteristics  using  data  from 
Aerospace  Corporation’s  Broadband  Array  Spectrograph  System  (BASS).  BASS  was  a 
non-imaging  spectrograph  that  operated  in  the  MWIR  and  LWIR  regions  of  the 
electromagnetic  spectrum.  The  Spatially  Enhanced  Broadband  Array  Spectrograph 
System  (SEBASS)  was  developed  as  a  follow-on  to  BASS.  SEBASS  is  an  imaging 
hyperspectral  sensor  that  was  "...  designed  to  explore  the  utility  of  infrared  hyperspectral 
sensors  for  remotely  identifying  solids,  liquids,  gases,  and  chemical  vapors  in  the  2.5  to 
14.0  p,m  'chemical  fingerprint'  spectral  region."  (Hackwell,  et  al,  1996) 

A.  THERMAL  INFRARED  MAPPING  SCANNER  (TIMS) 

The  Thermal  Infrared  Mapping  Scanner  (TIMS)  was  developed  for  remote 
sensing  of  non -renewable  resources.  TIMS  is  a  six  band,  dispersive  grating  scanner 
developed  by  Daedalus  Enterprises  for  the  Jet  Propulsion  Laboratory  (Richards,  1993). 
The  spectral  resolution  of  TIMS  was  designed  to  measure  the  depth  and  position  of  the 
reststrahlen  band  and  other  silicate  spectral  features  useful  in  identifying  rock  and  soil 
types.  The  spectral  bandpass  characteristics  of  TIMS  are  listed  in  Table  3.1. 
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Band  Number 

Spectral  Bands  (pm) 

1 

8.2  -  8.6 

2 

8.6  -  9.0 

3 

9.0  -  9.4 

4 

9.5-10.2 

5 

10.2-11.2 

6 

11.2-12.2 

Table  3.1.  TIMS  Spectral  Bands.  From  Richards,  1993. 


Data  from  TIMS  have  demonstrated  the  ability  of  thermal  spectral  data  to  conduct 
terrain  categorization  and  materials  classification  (Richards,  1993).  In  order  to  remove 
the  temperature  dependence  of  LWIR  emitted  radiation  and  focus  on  the  emissivity 
characteristics  of  the  materials  within  the  TIMS  field-of-view,  a  decorrelation  stretch  is 
often  applied  to  the  TIMS  data.  A  decorrelation  stretch  is  based  on  the  principal 
components  transformation.  The  processing  technique  “...causes  the  spectral  differences 
between  surface  units  to  be  displayed  as  color  differences,  while  most  of  the  temperature 
variation  is  displayed  as  intensity  differences.”  (Kahle  and  Goetz,  1983) 

An  analysis  of  TIMS  data  collected  over  the  Cuprite  Mining  District  of  Nevada  by 
Kahle  and  Goetz  (1983)  is  an  excellent  example  of  material  classification  using  TIMS 
data.  Figure  3.1  is  a  three-band  color  composite  of  the  Cuprite  Mining  District  from 
Kahle  and  Goetz.  The  discrimination  between  materials  is  based  primarily  on  the  free- 
silica  content  of  the  ground  materials  and  its  effect  on  their  emissive  spectra.  Silicified 
rocks  appear  as  bright  orange.  Clay-rich  minerals  are  magenta  and  carbonate  rocks  have 
a  green-blue  tone. 

B.  THE  BROADBAND  ARRAY  SPECTROGRAPH  SYSTEM  (BASS) 

BASS  is  a  dispersive  hyperspectral  instmment  capable  of  producing  a  complete 
1 16-band  spectrum  from  2.5-14.0  pm  every  5  milliseconds  (Bongiovi,  Hackwell,  and 
Hayhurst,  1995).  The  data  generated  by  BASS  demonstrated  the  utility  of  infrared 
spectral  data  and  prompted  the  development  of  an  imaging  instrument  based  on  the  same 
design. 
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1.  Instrument  Design 


BASS  employs  two  spectrometers  to  cover  the  MWIR  and  LWIR  regions  of  the 
electromagnetic  spectrum  (Bongiovi,  Hackwell,  and  Hayhurst,  1995).  The  MWIR 
channel  records  radiation  from  2.5  to  6.5  pm.  The  LWIR  channel  records  radiation  from 
6.5  to  14.0  pm.  Each  channel  employs  a  spherical-faced  prism  and  a  mirror  to  disperse 
and  reimage  the  incoming  radiation  onto  separate  58-element  linear  array 
photoconductors.  The  arrays  are  doped  silicon  ’blocked  impurity  band’  photoconductors 
which  have  a  high  quantum  efficiency  and  low  noise  when  operated  at  temperatures 
below  lOK.  The  optical  layout  of  BASS  is  presented  in  Figure  3.2.  Signals  from  the 
detectors  are  digitized  and  stored  by  a  Macintosh  n  personal  computer. 

BASS  has  been  flown  aboard  a  DH-6  "Twin  Otter"  aircraft.  The  BASS  cryostat 
is  mounted  on  a  platform  that  can  be  steered  in  the  roll  direction  of  the  aircraft.  A 
telescope  with  a  focal  length  of  one  meter  focuses  the  incoming  radiation  on  BASS’s  field 
optics.  The  field-of-view  of  the  telescope  is  one  milliradian.  From  operating  altitudes  of 
700  to  3000  ft,  the  ground  sample  distance  (GSD)  of  BASS  is  0.7  to  3.0  ft.  Three 
cameras  are  mounted  on  the  “Twin  Otter”  to  assist  the  operator  in  steering  BASS  and  to 
aid  in  post-mission  analysis  of  collected  spectra.  One  camera  is  a  wide  field-of-view 
color  camera  pointed  forward  at  a  45°  angle.  A  MWIR  camera  is  located  on  the  bottom 
the  platform  and  provides  a  8°  field-of-view.  A  black  and  white  camera  with  a  5°  field- 
of-view  is  bore  sighted  with  BASS  through  a  beam  splitter  in  the  external  telescope. 
Figure  3.3  is  a  schematic  of  the  BASS  "Twin  Otter"  configuration. 

2.  SMO  Utility 

An  analysis  of  BASS  hyperspectral  data  by  Bongiovi,  et  al  (1995)  demonstrated 
the  potential  of  LWIR  data  to  locate  and  identify  military  targets.  BASS  was  used  to 
collect  overhead  LWIR  hyperspectral  data  on  several  military  vehicles  in  October  of 
1994.  Data  from  the  BASS  collection  was  used  to  derive  spectral  indices  for  the  military 
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vehicles.  These  indices  were  then  applied  to  other  data  sets  to  determine  BASS’s  ability 
to  detect  target  vehicles  within  a  scene. 

Figure  3.4  shows  a  row  of  military  vehicles  that  were  used  to  derive  the  spectral 
templates.  BASS  data  collected  over  this  scene  was  converted  to  surface  brightness 
temperatures  after  atmospheric  compensation.  The  converted  data  was  then  used  to 
derive  four  spectral  templates:  sand,  limestone,  vehicle  1,  and  vehicle  2.  The  two  vehicle 
templates  were  representative  of  different  types  of  paint  observed  on  the  vehicles.  These 
templates  are  shown  in  Figure  3.5.  The  application  of  the  spectral  templates  to  new  data 
sets  results  in  response  plots  similar  to  Figures  3.6  and  3.7.  Figure  3.6  depicts  the 
response  plots  resulting  from  application  of  the  spectral  templates  to  the  initial  training 
scene.  The  arrows  indicate  a  response  to  the  presence  of  military  vehicles.  The 
application  of  the  spectral  templates  to  other  data  sets  demonstrated  an  ability  to  locate 
target  vehicles  even  when  partially  concealed  by  vegetation.  Figure  3.7  is  the  response 
function  of  a  tank  concealed  in  trees. 

C.  THE  SPATIALLY  ENHANCED  BROADBAND  ARRAY 
SPECTROGRAPH  SYSTEM 

The  success  of  BASS  led  to  the  development  of  an  imaging  infrared  instrument 
based  on  BASS’s  dispersive,  two-channel  design.  Aerospace  developed  the  follow-on 
instrument,  SEBASS,  with  internal  funds  as  a  phenomenology  research  tool  to  further 
LWIR  hyperspectral  remote  sensing  research. 

1.  Instrument  Design 

The  SEBASS  optical  layout  is  similar  to  that  of  BASS  and  is  depicted  in  Figure 
3.8.  Radiation  entering  SEBASS  is  split  into  two  spectrographs:  one  for  MWIR  (2.9-5.2 
p.m)  and  one  for  LWIR  (7.8-13.4  |j,m)  (Hackwell,  et  al,  1996).  Each  spectrograph  focuses 
a  dispersed  image  of  the  slit  onto  a  two-dimensional  array.  The  arrays  are  128  by  128 
element  impurity  band  conduction  arrays.  One  dimension  of  the  array  is  spectral  and  the 
other  is  spatial.  The  second  spatial  dimension  required  for  imaging  is  provided  by  the 
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movement  of  the  sensor.  Figure  3.9  displays  this  ‘pushbroom’  imaging  concept.  The 
figure  also  shows  the  spectral  and  spatial  dimensions  of  SEBASS  data. 

The  spectral  resolution  of  SEBASS  varies  across  the  array  due  to  the 
characteristics  of  longwave  infrared  radiation  and  the  dispersion  characteristics  of  the 
prisms.  Figure  3.10  is  a  graph  of  the  spectral  resolution  of  the  LWIR  channel.  SEBASS 
can  achieve  a  resolution  as  fine  as  .035  |im  at  the  long  wavelength  edge  of  the  array. 

SEBASS’s  data  system  collects,  digitizes  and  stores  signals  from  the  two  focal 
planes.  The  focal  planes  operate  at  data  rates  of  up  to  100  frames  per  second  and  are 
digitized  by  14-bit  analog-to-digital  converters.  Data  from  each  pass  is  stored  in  128 
megabytes  of  cache  RAM  until  after  a  pass  is  completed.  It  is  then  dumped  to  a  hard  disk 
for  later  processing.  The  system  is  built  around  a  Sun  Sparc  20  host  computer.  SEBASS 
can  be  mounted  aboard  the  Twin  Otter  aircraft  in  a  flight  configuration  similar  to  that 
used  during  BASS  flight  experiments.  SEBASS  can  employ  two  external  telescopes 
which  provide  either  a  1  milliradian  or  3  milliradian  field-of-view. 

2.  Data  Format  and  Use 

Prior  to  processing,  raw  SEBASS  sensor  data  must  be  spectrally  and 
radiometrically  calibrated.  Spectral  calibration  is  the  process  of  obtaining  an  accurate 
map  of  the  wavelength  band  collected  by  each  pixel  in  the  array.  Radiometric  calibration 
is  the  process  of  converting  instrument  digital  data  to  absolute  radiance  values. 

Data  produced  by  SEBASS  are  stored  according  to  the  following  convention: 

DN  {spectral,  spatial,  frame) 

DN  =  digital  number 

jpectra/  =  wavelength  dimension  (see  Figure  3.9) 
spatial  =  slit  position  dimension  (see  Figure  3.9) 
frame  =  ground  path  dimension  (see  Figure  3.9) 
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Spectral  Calibration 


a. 


The  wavelength  calibration  process  for  SEBASS’s  LWIR  channel  is 
performed  using  a  set  of  polymer  films  with  known  absorption  feature  (Bongiovi,  1996). 
The  transmittance  spectmm  of  the  filter  is  presented  as  Figure  3.11.  A  wavelength  map 
for  the  array  is  created  by  locating  the  filter’s  prominent  absorption  features  as  they  are 
imaged  on  the  focal  plane. 

The  polymer  film  is  placed  over  the  SEB  ASS  aperture  as  the  instmment 
images  a  blackbody  at  a  known  temperature.  The  blackbody  is  then  imaged  without  the 
filter.  The  ratio  of  these  two  images  is  the  transmittance  spectmm  for  the  filter.  Figure 
3. 12  is  an  example  of  the  spectmm  of  the  filter  as  measured  by  SEBASS.  The  location  of 
the  absorption  features  on  the  array  are  used  to  generate  a  polynomial  function  which 
describes  the  wavelength  for  each  pixel.  The  algorithm  which  generates  the  polynomial 
function  relies  on  the  fact  that  the  dispersion  of  SEBASS's  NaCl  prism  is  a  function  of 
the  square  of  wavelength  and  that  the  spatial  curvature  of  the  slit  is  a  second  order 
function.  The  final  fit  for  Aerospace's  LWIR  channel  wavelength  calibration  is  presented 
in  Equation  3.1. 

Equation  3.1  'k^{x,y)  =  aQ{y)  +  a^{y)x  +  a^{y)x^ 

flo  (y)  =  57.0758  -  0.0216030y  +  0.000100471y^ 
a,  (y)  =  1.03258  + 1.88537  x  10'^  y  +  9.2521 1x10"*  y^ 
a^{y)  =  -0.000100671 
X  =  spectral  dimension  (see  Figure  3.8) 
y  =  slit  position  dimension  (see  Figure  3.8) 

These  wavelength  calibration  equations  generate  a  wavelength  map  for  the 
LWIR  channel.  Due  to  the  curvature  of  SEBASS's  slit  and  dispersion  through  the  prism, 
the  center  wavelength  is  not  constant  across  the  array.  The  final  step  in  spectral 
calibration  is  the  removal  of  this  spectral  curvature.  This  is  accomplished  through  a 
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linear  interpolation  in  which  the  wavelength  map  of  SEBASS  is  interpolated  to  match  the 
center  column  of  the  array. 

h.  Radiometric  Calibration 

The  process  for  converting  instrument  digital  number  to  actual  radiance 
values  is  accomplished  through  the  use  of  on-board  blackbodies.  During  a  flight, 
SEBASS  collects  data  from  two  blackbodies  which  are  at  different  temperatures.  The 
calibration  sources  radiate  a  known  amount  of  energy  and  are  used  to  determine 
parameters  for  the  radiometric  calibration  of  SEBASS  data. 

The  focal  plane  array  employed  in  SEBASS  exhibited  a  non-linear 
response  function.  To  account  for  this  effect,  B.  Johnson  (1996)  developed  a  non-linear 
calibration  formula  which  is  presented  as  Equation  3.2.  The  calibration  files  collected 
from  the  blackbodies  are  used  to  generate  a  series  of  coefficients  {BG,a,x)  which 
convert  the  instrument’s  digital  values  into  measured  radiance ,  .  The  quantities  a  and 

X  are  fixed  parameters  for  the  SEBASS  instmment  during  a  collection  exercise.  The 
BG  parameter  represents  the  background  system  noise  which  may  vary  from  flight  to 
flight. 


Equation  3.2  ^Vcoadds  =  +  BG)  +x(L„  -H  BGf 

L„  =  measured  radiance  / cm^\m  sr\ 

BG=  background  radiance  [p,W/cm^pm  5r] 

c,T  =  gain  coefficients 

DN  =  instrument  digital  number 

Coadds  —  number  of  co-added  measurements 

The  initial  configuration  of  SEBASS  exhibited  a  small  degree  of  thermal 
drift.  To  account  for  this  effect,  the  background  radiance  term,  BG ,  for  a  particular  data 
set  is  adjusted  by  interpolating  between  the  two  closest  calibration  measurements.  This 
interpolation  is  accomplished  through  Equation  3.3. 


31 


Equation  3.3 


5G(0  =  BG,+(5G2-5G, 


BG  =  background  radiance  /cm  ^  |Lim  ^r] 

t=time  [sec] 

a  =0.038 
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IV.  EXERCISE  WESTERN  RAINBOW 


The  MASINT  Hyperspectral  Study  (MHS)  was  initiated  in  1993  to  address  the 
utility  of  hyperspectral  imaging  systems.  The  MHS  program  conducted  a  series  of 
experiments  focusing  on  the  ability  of  hyperspectral  imaging  technology  to  address  the 
following  intelligence  areas: 

•  Support  to  Military  Operations  (SMO) 

•  Nuclear  Proliferation  (MHS:  Interim  Progress  Report,  1995) 

The  MHS  program  adopted  a  standardized  experimental  method  for  its 
hyperspectral  collection  exercises.  The  MHS  exercises  combine  realistic,  military 
scenarios  with  a  succession  of  laboratory  measurements,  ground  field  measurements,  and 
airborne  collections. 

During  1994  and  1995,  the  MHS  program  conducted  two  experiments  to  study  the 
SMO  capabilities  of  hyperspectral  imagery.  In  September  1994,  Exercise  SOUTHERN 
RAINBOW  was  conducted  at  the  Chicken  Little  facility  at  Eglin  Air  Force  Base.  A 
follow-on  exercise,  WESTERN  RAINBOW,  was  conducted  in  October  1995  at  the  Yuma 
Proving  Ground  in  Yuma,  Arizona  (MHS:  Interim  Progress  Report,  1995).  The  primary 
objective  of  WESTERN  RAINBOW  was  "...  to  measure  the  multi/hyperspectral  and 
radar  signatures  of  ground  threat  vehicles  and  backgrounds  under  various  conditions  of 
concealment,  countermeasures,  and  lighting."  (Operations  Order,  1995) 

A.  EXERCISE  SCENARIO 

In  order  to  provide  a  realistic  environment  for  target  detection  and  classification,  a 
scripted  scenario  was  used  to  guide  the  selection,  placement,  and  camouflage  of  military 
targets  used  in  WESTERN  RAINBOW.  Personnel  from  the  Joint  Camouflage, 
Concealment,  and  Deception  (JCCD)  group  of  the  Office  of  the  Secretary  of  Defense 
drafted  the  exercise  scenario  and  directed  the  collection  site  preparations.  The  scenario 
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defined  the  placement  of  tactical  military  vehicles  as  well  as  their  camouflage  and 
concealment  schemes.  The  JCCD  efforts  ensured  that  the  deception  schemes  used 
current  camouflage,  concealment,  and  deception  (CCD)  equipment  and  techniques. 

The  scripted  scenario  for  WESTERN  RAINBOW  described  a  decimated  orange 
force  defending  critical  mobile  targets  (CMTs)  from  a  southward  moving  blue  force 
(WESTERN  RAINBOW:  CCD  Report,  1996).  Two  sites  with  multiple  levels  of  CCD 
were  prepared.  A  tank  and  artillery  site  (referred  to  as  the  JCCD  site)  was  established  to 
represent  a  forward  defense  position.  A  second  site  (referred  to  as  the  Arroyo  site) 
represented  the  orange  forces’  CMTs. 

The  JCCD  site  consisted  of  four  Soviet  tanks  and  a  mobile  artillery  piece.  In  the 
scenario,  this  site  was  established  with  the  mission  of  defending  the  critical  targets  from 
the  southward  moving  blue  forces.  Nine  decoys  of  varying  fidelity  were  placed  just  north 
of  the  five  defending  vehicles.  These  decoys  were  designed  to  have  a  greater  signature 
visibility  than  the  real  targets  in  order  to  reduce  the  chance  of  effective  aerial 
reconnaissance  or  targeting  by  the  blue  forces. 

The  Arroyo  site  represented  the  orange  forces  theater  ballistic  missile  (TBM) 
assets  and  their  associated  close-in  protection.  Three  Soviet  tanks  and  a  Soviet  armored 
troop  carrier  were  deployed  to  protect  the  two  TBMs.  The  tactical  vehicles  were  placed 
in  defensive  positions  around  the  TBMs. 

B.  COLLECTION  SITES 

The  SMO  portion  of  WESTERN  RAINBOW  consisted  of  three  collection  sites. 
Each  collection  site  served  a  specific  purpose  within  the  overall  collection  plan.  The 
Malpai  target  site  served  as  an  initial  target  and  background  signature  collection  site. 

This  site  was  utilized  for  spectral  ground  truth  collection  and  signature  characterization. 
The  JCCD  tank  and  artillery  site  and  Arroyo  TBM  site  were  established  for  realistic  target 
detection  analysis. 
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1.  Malpai  Site 


,  The  Malpai  site  consisted  of  an  array  of  threat  vehicles,  camouflage  schemes, 
decoys,  and  calibration  panels.  The  layout  of  the  targets  within  the  Malpai  site  is 
depicted  in  Figure  4.1  (Boisjolie,  1996).  The  array  was  situated  on  a  background  of 
‘desert  pavement’  which  consisted  of  small  pebbles  coated  with  a  desert  varnish. 
Vegetation  was  sparse  and  located  in  two  washes  that  ran  alongside  the  collection  site. 
The  site  was  established  on  19  October  and  maintained  for  five  days.  The  array  was 
modified  periodically  by  the  addition  of  camouflage  and  the  removal  or  addition  of 
various  tactical  vehicles.  These  changes  were  documented  for  each  collection  mission. 

An  array  of  five  Soviet  T-72  tanks  was  centered  within  the  target  site.  The  tanks 
were  oriented  in  various  directions  to  allow  for  multiple  sun-angle  geometries.  Two  US 
tankers  were  placed  to  the  west  of  the  T-72  array.  In  order  to  provide  data  for  thermal 
inertia  and  temperature  analysis,  one  tank  was  filled  with  water  while  the  other  was  left 
empty.  Two  sets  of  reflective  calibration  panels  were  placed  within  the  array  for  the 
calibration  of  visible  and  NIR  sensors.  In  the  north  part  of  the  site,  a  series  of  JCCD 
camouflage  treatments  and  decoys  were  placed.  These  included  US  and  Soviet  color 
panels,  a  special  IR  applique,  prototype  camouflage  materials,  and  helicopter  decoys.  At 
the  south  end  of  the  site,  a  series  of  threat  vehicles  were  placed  in  the  open.  A  general  list 
of  targets  within  the  Malpai  site  is  presented  in  Table  4. 1 . 
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Target 

Bumper  # 

T-72 

199 

T-72 

934 

T-72 

922 

T-72 

910 

T-72 

933 

2S3 

5521 

T-80  surrogate 

001 

T-80  surrogate 

002 

LAV  (green) 

5217 

LAV  (desert  camouflage) 

82 

M  series  5000  gal  tanker  (empty) 

Cl  38 

M  series  5000  gal  tanker  (filled  with  water) 
George 

Calibration  panels 

Calibration  panels 

Calibration  pool  (filled  with  water) 

Apache  helicopter  decoy  (2 1/2  dimension  mural) 
OH-58  helicopter  inflatable  decoy 

PRC  net 

Bomb  crater  decoy 

DOD  Less  arctic  net 

US  panel  of  colors 

Soviet  panel  of  colors 

Prototype  fox  hole  camouflage 

Special  IR  applique 

Prototype  ULCAN 

Prototype  tank  mat 

C243 

Table  4. 1 .  Malpai  Target  Descriptions.  After  Boisjolie,  1996. 


2.  JCCD  Tank  and  Artillery  Site 

The  JCCD  tank  and  artillery  site  employed  multiple  threat  vehicles  and  decoys 
with  varying  levels  of  CCD.  The  site  consisted  of  three  rows  of  targets  which  were 
placed  within  washes  running  from  the  southwest  to  the  northeast.  The  layout  of  the 
JCCD  site  is  depicted  in  Figure  4.2  (Boisjolie,  1996).  The  washes  contained  sparse 
vegetation  consisting  of  small  bushes  and  trees.  The  JCCD  site  was  maintained  from  19 
October  through  24  October. 

Two  of  the  target  lines  contained  decoys  which  varied  in  quality  from  high 
fidelity,  specially  constmcted  mock-ups  to  field  expedient  decoys  constructed  from 
available  scraps  found  within  the  target  site.  Additional  deception  devices  such  as  comer 
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reflectors,  heating  pots,  and  fake  tank  tracks  were  used  to  enhance  the  decoys.  A  full 
description  of  each  decoy  is  included  in  Table  4.2. 


Tarqet 

Label  in  Figure  4.3 

Additional  Description 

M-19  aluminum  matting 
tank  decoy 

Decoy  1 ,  Line  1 

No  camouflage 

Sheet  metal  tank  decoy 

Decoy  2,  Line  1 

Partially  covered  by  LOSS  desert 
camouflage  netting 

M-19  aluminum  matting 
tank  decoy 

Decoy  3,  Line  1 

Covered  with  AFCENT  camouflage 
netting 

T-72  mural 

Decoy  4,  Line  1 

Less  army  woodland  camouflage 
netting  placed  beside  mural,  heated 
thermal  reflector 

M-1  turret  mural 

Decoy  5,  Line  1 

Less  army  woodland  camouflage 
netting  underneath,  fresh  cut 
vegetation 

2S3  sled  decoy 

Decoy  1 ,  Line  2 

Covered  with  AF  ULCAN  desert 
camouflage  netting,  propane 
heaters 

2S3  sled  decoy 

Decoy  2,  Line  2 

Covered  with  LCSS  army  woodland 
camouflage  netting 

Scrap  metal  tank  decoy 

Decoy  3,  Line  2 

Covered  with  West  German 
woodland  camouflage  netting, 
propane  heaters 

M-1  mural 

Decoy  4,  Line  2 

East  German  camouflage  netting 
placed  in  front 

Table  4.2.  JCQD  Decoy  Line  Descriptions.  After  Boisjolie,  1996. 


The  real  threat  vehicles  were  deployed  within  the  southern  wash.  These  vehicles 
were  camouflaged  and  placed  in  both  full  and  partial  defilade.  The  full  dress  for  each 
vehicle  is  outlined  in  Table  4.3. 


Target 

Bumper  # 

Label  in 
Figure  4.3 

Additional  Description 

2S3  mobile  artillery 

5521 

Artillery 

Camouflage  with  cut  live  and  dead 
vegetation 

T-72 

910 

Tank  4 

Camouflage  with  live  &  dead 
vegetation  and  dirt 

T-72 

199 

Tank  3 

Full  defilade,  covered  with  sandbags 
and  Less  desert  camouflage  netting 

T-72 

934 

Tank  2 

Camouflage  with  dead  vegetation  and 
soil 

T-72 

933 

Tank  1 

Partial  defilade,  UK  desert 
camouflage  netting,  UK  thermal 
blanket 

Table  4.3.  JCCD  Vehicle  Target  Line  Description.  After  Boisjolie,  1996. 


3.  Arroyo  Site  TBM  Site 

The  Arroyo  site  was  a  simulation  of  a  TBM  site  and  its  associated  defensive 
weapons.  The  site  had  a  moderate  amount  of  vegetation.  The  Arroyo  site  was 
established  on  21  October  and  maintained  until  23  October.  The  layout  for  the  Arroyo 
site  is  presented  in  Figure  4.3  (Boisjolie,  1996).  The  targets  included  two  TBMs,  three 
defensive  tanks,  and  a  personnel  carrier.  All  of  the  targets  within  the  Arroyo  utilized 
some  kind  of  camouflage  treatment.  The  complete  description  of  the  CCD  treatments  for 
each  target  are  listed  in  Table  4.4. 


Target 

Bumper  # 

Label  in 
Figure  4.4 

Additional  Description 

T-80  surrogate 

002 

Tank  1 

AF  woodland  ULCAN  camouflage  netting 

T-80  surrogate 

001 

Tank  2 

WSMR  just-in-time  netting 

T-72 

922 

Tank  3 

Partial  defilade,  Isreali  ULCAN  netting 

MT-LB 

PM  1-904 

Tank  4 

Partial  defilade,  SSDC  space  frame 

TBM  1 

TBM1 

TBM  1 

Canvas,  LCSS  tri-color  desert 

camouflage  netting 

TBM  2 

TBM  2 

TBM  2 

Woodland  Eglin  ULCAN  netting 

Table  4.4.  Arroyo  Target  Description.  After  Boisjolie,  1996. 
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C.  GROUND  COLLECTION  EQUIPMENT 


The  ground  sensors  employed  in  WESTERN  RAINBOW  provided  an  extensive 
data  set  of  background  and  target  spectra.  The  complete  ground  sensor  suite  is  listed  in 
Table  4.5.  As  indicated  in  the  table,  the  electromagnetic  spectrum  from  visible  to  LWIR 
was  covered.  Instrumentation  capable  of  collecting  MWER  and  LWIR  data  during  the 
exercise  included  the  Johns  Hopkins  University  (JHU)  imcro  fourier  transform  imaging 
radiometer  (uFTIR)  and  the  Environmental  Institute  of  Michigan’s  (ERJM)  Bomem  MB- 
100  FTIR  (Operations  Order,  1995).  In  addition,  JHU  performed  laboratory  MWIR  and 
LWIR  measurements  of  material  samples  with  a  Nicolet  FTIR  spectrometer  and  a 
Beckman  grating  spectrometer. 


Instument 

Spectral  Range  (urn) 

Spectral  Resolution  (cm  -1 ) 

GER  Mark  IV 

0.35-2.5  (875  bands) 

0.005 

JHU  micro-FTIR 

3-5.5, 7.5-14 

6 

Bomem  MB-100  FTIR 

3-15  (128  bands) 

1-4 

Nicolet 

2-15 

4 

Beckman 

.4-2.5 

2nm  visible,  10  nm  SWIR 

Table  4.5.  WESTERN  RAINBOW  Ground  Sensors.  After  Operations  Order,  1995. 


D.  AIRBORNE  COLLECTION  EQUIPMENT 

The  airborne  collection  assets  involved  in  WESTERN  RAINBOW  covered  the 
electromagnetic  spectrum  from  visible  to  LWIR  and  included  both  multispectral  and 
hyperspectral  sensors.  The  sensors  are  listed  in  Table  4.6.  The  assets  capable  of  LWIR 
remote  sensing  included  Aerospace  Corporation’s  SEBASS  and  the  EG&G's  Deadalus 
multispectral  scanner  (Daedalus  MSS)  (Operations  Order,  1995).  The  Daedalus  MSS  is  a 
multispectral  scanner  with  12  bands.  One  channel  can  be  configured  for  MWIR  imaging 
and  one  channel  can  be  configured  for  LWIR  imaging. 
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Instrument 

Spectral  Range(|i,m) 

Number  of  Spectral 
Bands 

Ground  Sample 
Distance  (feet) 

HYDICE 

.4-2.5 

210 

2.5-10 

Daedalus  MSS 

.4-12 

12  (8  vis-NIR,  2 
SWIR,  1  MWIR,  1 

0.5 

SEBASS 

7-14 

LWIR) 

128 

1-3 

Table  4.6.  WESTERN  RAINBOW  Airborne  Sensors.  After  Operations  Order,  1995. 
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V.  EXPERIMENT  OBSERVATIONS 


The  WESTERN  RAINBOW  experiment  was  conducted  from  18-24  October, 
1995  at  the  United  States  Army  Proving  Ground  in  Yuma,  Arizona.  Airborne  SEBASS 
data  collection  flights  over  the  SMO  target  areas  were  conducted  from  20  -  24  October. 
The  site  preparations,  ground  truth  observations,  and  airborne  sensor  operations  were 
coordinated  to  ensure  that  the  target  areas  were  well  documented  for  post-exercise 
analysis.  This  chapter  presents  the  author’s  observations  of  the  exercise.  The  chapter 
begins  with  a  description  of  the  Yuma,  Arizona  collection  environment.  The  second 
section  describes  the  performance  of  the  SEBASS  sensor.  The  final  section  describes  the 
hyperspectral  data  collected  by  SEBASS  and  the  procedures  used  to  prepare  it  for 
analysis. 

A.  YUMA  PROVING  GROUND  OBSERVATIONS 

The  Yuma  Proving  Ground  was  an  excellent  environment  for  the  exercise.  The 
area  has  predictable  weather  conditions  and  low  humidity.  The  airspace  over  the  site  is 
restricted.  This  allowed  sensor  aircraft  to  maneuver  in  the  absence  of  commercial  traffic. 
The  average  elevation  was  about  850  feet. 

1.  Collection  Sites 

The  isolation  of  the  exercise  area  provided  an  environment  in  which  the 
collection  sites  could  be  established  and  camouflaged  without  interference.  The  desert 
background  consisted  of  sand,  rock  fragments,  and  small  brush.  The  rock  fragments  were 
darkened  with  a  distinctive  desert  varnish.  Figure  5.1  is  a  picture  of  the  JCCD  area 
displaying  the  general  background  of  the  Yuma  Proving  Ground. 

The  collection  sites  were  prepared  by  the  JCCD  personnel  according  to  the 
scenario  described  in  Chapter  IV,  Exercise  WESTERN  RAINBOW.  The  camouflage  and 
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concealment  preparations  were  realistic  and  simulated  actual  operating  conditions.  The 
site  preparations  and  changes  were  well  documented  and  photographed. 

Ground  traffic  in  the  collection  sites  was  kept  to  a  minimum  with  the  exception  of 
the  Malpai  site  which  was  used  extensively  for  ground  sensor  collections.  In  addition, 
traffic  in  the  JCCD  site  was  limited  to  ground  truth  personnel  in  order  to  provide  a  “blind 
test”  of  sensor  performance. 

2.  Weather  Conditions 

The  weather  conditions  were  consistent  throughout  the  exercise  and  were 
characterized  by  hot,  dry  days.  The  surface  temperature  during  the  exercise  ranged  from 
lows  of  about  20°  C  at  night  to  as  high  as  54°  C  during  the  day.  Wind  conditions  were 
generally  good  throughout  the  exercise.  Only  one  day  of  operations  had  to  be  canceled 
because  of  low  visibility  caused  by  wind.  Meteorological  data  such  as  wind  speed, 
surface  temperature,  air  temperature,  and  humidity  were  collected  by  ground  truth 
personnel. 

B.  SEBASS  PERFORMANCE 

The  Aerospace  SEBASS  instrument  was  flown  aboard  a  Twin  Otter  aircraft  in  a 
configuration  similar  to  Figure  3.3.  The  author  observed  both  collection  and  post-mission 
processing  of  SEBASS  data.  The  SEBASS  flights  were  flown  at  0600,  1000,  and  1400 
local  time  on  various  days  during  the  exercise.  Each  flight  consisted  of  several  runs  over 
the  various  collection  sites.  Coordination  between  the  aircraft  and  on-going  ground 
operations  was  accomplished  through  daily  briefings  and  radio  contact. 

1.  Flight  Operations 

The  SEBASS  flights  over  the  Arroyo,  JCCD,  and  Malpai  collection  sites  were 
performed  from  an  altitude  of  2,500  feet.  This  altitude  produced  images  with  a  ground 
pixel  size  of  2.5  feet.  The  swath  width  of  each  data  set  was  approximately  320  feet. 
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The  small  swath  width  of  SEBASS  contributed  to  difficulties  in  getting  a  good 
image  of  all  of  the  targets  within  a  collection  site.  SEBASS  was  installed  aboard  the 
aircraft  without  roll  stabilization.  As  a  result,  wind  gusts  caused  targets  to  shift  in  and  out 
of  the  instrument  s  field-of- view  as  the  aircraft  rolled.  In  addition,  the  lack  of  ground 
markers  and  on-board  positioning  equipment  made  it  very  difficult  for  the  pilot  to  fly 
directly  over  a  particular  set  of  targets  and  keep  them  within  SEBASS’s  small  field-of- 
view.  The  additional  cameras  aboard  the  aircraft  and  the  sensor  operator’s  ability  to  steer 
the  instrument  partially  alleviated  this  problem.  White  panel  markers  and  ground-to-air 
voice  communications  were  used  on  later  flights  to  help  steer  the  pilot  over  a  particular 
target  area. 


2.  Sensor  Performance 

Exercise  WESTEN  RAINBOW  was  the  first  operational  use  of  SEBASS.  The 
data  collected  served  as  the  first  image  data  recorded  by  the  instrument  outside  of  the 
laboratory.  The  instrument  performed  extremely  well  throughout  the  exercise  and 
maintained  a  noise-equivalent  spectral  radiance  (NESR)  1.2  x  10“®  W/cm^pm  sr . 

Several  artifacts  were  noticed  in  the  collected  data.  Intensely  bright  objects 
created  a  ghost  image  of  themselves  approximately  ten  lines  ahead  of  their  actual 
location.  Figure  5.2  shows  an  example  of  this  artifact.  Each  T-72  in  the  figure  has  a 
ghost  image  offset  by  several  scan  lines.  The  problem  was  caused  by  the  reflection  and 
diffraction  of  incoming  radiation  as  it  interacted  with  the  edge  of  the  imaging  slit.  The 
ultimate  effect  of  the  ghost  problem  was  a  mixing  of  scene  radiance  from  each  line  with 
the  line  ten  scans  ahead  of  it.  This  artifact  did  not  have  a  significant  effect  on  the  results 
of  this  study.  A  second  problem  identified  during  the  exercise  was  the  thermal  drift 
mentioned  in  Chapter  m.  Section  C.  The  thermal  drift  was  initially  detected  and 
characterized  using  data  collected  during  WESTERN  RAINBOW.  Finally,  the  gain  of 
one  of  the  four  amplifiers  in  SEBASS  varied  from  flight  to  flight.  Although  these 
problems  have  since  been  corrected  in  the  instrument,  they  affect  all  of  the  data  collected 
during  WESTERN  RAINBOW  to  some  degree. 
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c. 


GROUND  TRUTH  OBSERVATIONS 


The  collection  of  ground  truth  target  and  background  spectra  was  accomplished 
by  ERIM  and  John  Hopkins  University  personnel  using  the  equipment  listed  in  Table  4.5. 
The  ERIM  collection  operations  were  performed  using  an  elevated  platform.  The  John 
Hopkins  University  collections  were  performed  using  the  micro-FTIR,  Nicolet,  and 
Beckman  spectrometers. 

The  John  Hopkins  University  micro-FTIR  employed  in  the  field  recorded 
measurements  in  radiance  values.  These  measurements  were  converted  to  apparent 
emissivity  for  spectral  analysis  of  the  SEB  ASS  data.  The  conversion  process  began  with 
the  estimation  of  the  surface  temperature  of  each  sample.  The  sample  surface 
temperature  was  estimated  by  calculating  the  brightness  temperature  within  each 
wavelength  band.  Brightness  temperature  was  calculated  by  inverting  Equation  2.1  and 
setting  the  surface  temperature  equal  to  the  maximum  brightness  temperature  measured 
within  the  bandwidth  of  the  instrument.  This  estimation  is  based  on  the  assumption  that 
the  emissivity  maximum  occurs  somewhere  between  8.0  and  13.0  p,m  and  is  close  to  1.0. 
The  apparent  emissivity  of  the  material  was  calculated  through  application  of  Equation 
2.2. 

Ground  truth  spectra  were  available  on  several  targets  and  camouflage  schemes. 
The  available  spectra  and  their  source  are  listed  in  Table  5.1  and  depicted  graphically  in 
Figures  5.3  and  5.4.  For  those  materials  for  which  ground  truth  spectra  were  not 
available,  the  author  sampled  spectra  directly  from  SEB  ASS  images  of  the  Malpai 
training  scene.  The  Malpai-derived  spectra  are  listed  in  Table  5.2.  Figure  5.5  depicts 
some  examples  of  spectra  extracted  from  SEBASS  data. 
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Target 

Source 

Designator 

Additional  Description 

Desert  Varnish 

JHU 

MAL1048_JHU 

Sand  Beneath  Desert  Varnish 

JHU 

DUST1110_JHU 

Desert  Panorama 

JHU 

DSRT1606_JHU 

T-72  Tank 

JHU 

T1326_JHU 

T-80  Surrogate 

JHU 

MBT1152_JHU 

Woodland  Camouflage 

JHU 

WOOD1347_JHU 

British/German 

Less  Desert  Camouflage 

JHU 

LC1531_JHU 

Desert  Tricolor 

LCDD  Woodland  Camouflage 

JHU 

LC1322_JHU 

Woodland  Pattern 

Tan  Camouflage 

JHU 

TAN1305_JHU 

AF  Desert  Ulcan 

Table  5. 1 .  Ground  Truth  Spectra  from  WESTERN  RAINBOW. 


Target 

Source 

Desert  Varnish 

SEBASS  -  20oct95_1000/mal5.dat 

Disturbed  Desert  Varnish 

SEBASS  -  20oct95_1 000/mal5.dat 

Tank  Tracks  in  Desert  Varnish 

SEBASS  -  20oct95_1000/mal5.dat 

Vegetation 

SEBASS  -  20oct95_1000/mal5.dat 

Prototype  Foxhole  Cover 

SEBASS  -  20oct95_1000/mal5.dat 

IR  Applique 

SEBASS  -  20oct95_1000/mal5.dat 

Table  5.2.  Target  Spectra  Derived  from  SEBASS  Data  Collected  Over  Malpai  Site. 


D.  HYPERSPECTRAL  OBSERVATIONS 

Data  from  WESTERN  RAINBOW  was  supplied  to  the  author  by  the  Office  of 
Research  and  Development,  Aerospace  Corporation,  and  John  Hopkins  University.  The 
amount  of  image  data  from  SEBASS  was  extensive.  Each  SEBASS  image  is 
approximately  140  megabytes.  As  many  as  14  images  were  collected  during  a  flight.  The 
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author  screened  the  SEB  ASS  data  set  and  selected  particular  collection  flights  for 
processing  and  analysis. 

1.  Hyperspectral  Image  Data 

Data  from  the  Malpai,  Arroyo,  and  JCCD  sites  were  analyzed  to  determine  if 
thermal  spectral  data  would  provide  an  ability  to  detect  military  targets  under  various 
degrees  of  concealment.  The  Malpai  site  was  used  as  an  initial  training  scene.  It  served 
as  a  source  of  target  spectra  and  as  a  tool  with  which  to  select  how  the  principal 
components  transformation  would  be  applied  to  later  data  sets.  The  JCCD  and  Arroyo 
scenes  were  then  analyzed  to  determine  if  the  various  targets  could  be  detected. 

The  selection  of  which  Arroyo  and  JCCD  scenes  to  analyze  was  driven  by  the 
requirement  to  maximize  the  possibility  of  an  image  set  containing  all  of  the  targets 
within  a  collection  site.  The  data  were  also  screened  to  determine  which  data  sets  were 
not  adversely  affected  by  the  amplifier  gain  problem.  The  images  listed  in  Table  5.3  were 
selected  for  analysis. 


Date 

Time 

Collection  Site 

Data  Set  Name 

20-Oct 

1000 

Malpai 

20oct95_1 000/mal5.dat 

23-Oct 

1000 

JCCD  -  Line  1 

23oct95_1 000/jccd3.dat 

23-Oct 

1000 

JCCD  -  Line  2 

23oct95_1 000/jccdl  .dat 

23-Oct 

1000 

JCCD  -  Vehicle  Target  Line 

23oct95_1 000/iccd4.dat 

23-Oct 

1000 

Arroyo 

23oct95_1 000/arroyo5.dat 

Table  5.3.  SEBASS  Data  Sets  Selected  for  Analysis. 


2.  Data  Preparation 

The  raw  SEBASS  data  were  processed  into  several  standard  data  formats.  The 
steps  involved  in  the  data  preparation  process  are  depicted  in  Figure  5.6. 
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a.  Spectral  and  Radiometric  Calibration 

The  spectral  and  radiometric  calibration  procedures  described  in  Chapter 
in,  Section  C  were  applied  to  the  raw  data  to  convert  the  data  into  appropriate  units.  The 
spectral  curvature  of  the  data  was  removed  in  this  step  through  linear  interpolation. 

The  effect  of  the  varying  gain  of  the  fourth  amplifier  was  mitigated 
through  a  procedure  proposed  by  B.  Johnson  (1996)  of  Aerospace  Corporation.  The 
output  of  the  fourth  amplifier  was  compared  to  that  of  its  adjacent  amplifier  outputs.  A 
multiplication  factor  was  then  derived  to  bring  the  output  of  the  errant  amplifier  down  to 
the  level  of  the  other  amplifiers.  This  procedure  was  based  on  the  assumption  that  the 
gain  of  the  fourth  amplifier  was  constant  across  the  focal  plane  array  and  the  taps  from 
the  array  were  staggered  enough  to  allow  noise  and  material  spectral  characteristics  to  be 
averaged  out.  The  SEBASS  data  was  corrected  for  gain  effects  prior  to  radiometric 
calibration. 

b.  Generation  of  a  Simulated,  FLIR  Image 

The  radiometrically  and  spectrally  calibrated  data  was  then  used  to 
generate  a  simulated,  single-band  FLIR.  This  image  was  generated  by  averaging  the 
radiance  values  measured  in  each  band.  It  served  as  the  standard  image  with  which  to 
compare  the  results  from  the  analysis  of  the  spectral  data 

c.  Atmospheric  Correction 

Atmospheric  correction  was  a  two  step  process.  The  plastic  ruler 
atmospheric  correction  procedure  was  used  to  determine  relative  atmospheric 
transmittance  and  upwelling  radiance.  These  results  were  then  scaled  using  Aerospace 
Corporation’s  ARTSPEC  atmospheric  model. 

The  plastic  ruler  algorithm  outlined  in  Chapter  11,  Section  B  was  modified 
slightly  to  account  for  the  effects  of  employing  the  algorithm  to  a  scene  composed  almost 
entirely  of  desert  sand.  The  plastic  ruler  algorithm  works  well  in  areas  with  large 
amounts  of  vegetation  because  vegetation  has  a  relatively  flat  emissive  spectrum  which  is 
close  to  1.0  throughout  the  LWIR  region.  Unfortunately,  the  desert  environment  in  Yuma 
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has  very  sparse  vegetation  and  most  of  the  pixels  within  each  image  were  desert  sand. 
The  emissive  spectrum  of  desert  sand  has  a  characteristic  reststrahlen  minimum  between 
8.0  and  10.0  |xm.  As  the  plastic  ruler  algorithm  attempts  to  fit  a  line  through  the  upper 
bound  of  a  scatter  plot  composed  almost  entirely  of  sand  measurements,  it  produces  an 
artificially  low  transmittance  in  the  8.0  to  10.0  |j,m  region.  To  solve  this  problem,  the 
author  screened  the  data  to  exclude  those  pixels  that  exhibited  a  strong  reststrahlen 
feature.  This  screening  was  performed  by  comparing  the  brightness  temperature  of  each 
pixel  between  9.45  and  10.0  |j.m  and  excluding  those  pixels  whose  difference  exceeded  a 
particular  threshold.  The  author  also  found  that  a  large  percentage  of  the  desert 
vegetation  exhibited  a  high  degree  of  senescence.  If  a  large  percentage  of  senescent 
vegetation  was  used  to  generate  the  atmospheric  parameters,  the  spectral  characteristics 
of  senescent  vegetation  described  in  Figure  2.7  could  also  adversely  affect  the  results.  To 
prevent  this  from  occurring,  pixels  exhibiting  a  large  degree  of  senescence  were  excluded 
by  comparing  the  brightness  temperature  between  10.65  and  1 1.5  jxm. 

The  plastic  ruler  atmospheric  compensation  algorithm  was  implemented  in 
DDL  routines  and  is  included  as  Appendix  B.  Figures  5.7  to  5.9  depict  the  results  of  the 
algorithm  for  SEB  ASS  data  collected  over  the  Malpai  site  at  1000  local  time  on  20 
October  1995.  Figure  5.7  depicts  a  region  of  the  scene  and  the  two  screening  stages  in 
which  the  pixels  representing  healthy  vegetation  are  selected.  The  image  sequence  shows 
the  pixels  which  satisfied  the  threshold  check  at  each  screening  stage.  Figure  5.8  depicts 
the  results  of  the  plastic  ruler  algorithm  if  no  screening  has  taken  place.  The  first  graph  is 
the  transmittance  spectrum  of  the  atmosphere.  The  second  graph  is  the  upwelling 
radiance  emitted  by  the  atmosphere  and  received  by  the  sensor.  The  effect  of  the  sand 
pixels  can  be  seen  as  a  dip  in  atmospheric  transmittance  between  9.0  and  10.0  pm. 

Figure  5.9  depicts  the  results  of  the  plastic  ruler  algorithm  after  the  input  data  has  been 
screened  to  include  only  healthy  vegetation. 

The  results  of  the  plastic  ruler  algorithm  were  then  compared  with  the 
results  of  Aerospace  Corporation’s  ARTSPEC  atmospheric  model.  The  plastic  ruler 
algorithm  uses  the  collected  data  to  generate  relative  atmospheric  parameters.  The 
ARTSPEC  atmospheric  model  was  used  to  convert  the  plastic  ruler  results  into  true 
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atmospheric  transmittance  and  upwelling  radiance.  The  comparison  was  performed  using 
a  series  of  routines  written  by  Burke  (1995).  The  routines  were  originally  written  for  the 
BASS  instrument  and  were  modified  by  the  author  to  work  with  SEBASS  data.  The 
routines  alter  the  input  parameters  of  the  ARTSPEC  model  until  the  model’s  output 
atmospheric  profile  matches  the  shape  of  the  results  from  the  plastic  ruler  algorithm. 
Figure  5.10  depicts  the  results  of  the  ARTSPEC  fitting  procedure  on  the  atmospheric 
parameters  generated  from  the  Malpai  data.  It  shows  the  ARTSPEC  model  results  as  a 
solid  curve  and  the  scaled  plastic  ruler  algorithm  results  as  a  dotted  curve. 

The  scaled  atmospheric  transmittance  and  upwelling  radiance  parameters 
were  then  applied  to  the  calibrated  data  sets  to  remove  the  effects  of  the  atmosphere.  This 
was  done  by  inverting  Equation  2.8. 

d.  Conversion  to  Apparent  Emissivity 

The  atmospherically  corrected  data  were  then  converted  into  apparent 
emissivity  using  the  same  procedure  that  was  applied  to  the  JHU  radiance  data  as 
described  in  Section  C  above.  Figure  5.1 1  compares  the  apparent  emissivity  of  desert 
varnish  extracted  from  the  Malpai  data  with  JHU  ground  truth  data  collected  from  the 
same  area.  The  two  ground  truth  spectra  represent  the  range  of  values  observed  in 
spectral  measurements  taken  from  several  different  desert  varnish  samples. 

3.  Hardware  and  Software  Configuration 

The  author  used  a  UNIX-based,  UltraSparc  Sun  workstation  for  data  analysis. 

The  three  hyperspectral  analysis  algorithms  were  implemented  in  the  ENVI  hyperspectral 
data  analysis  software  package.  ENVI  was  used  to  perform  the  principal  components 
analysis,  spectral  angle  mapper  analysis,  and  spectral  matched  filter  analysis. 
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VI.  DATA  ANALYSIS 


The  four  data  sets  identified  in  Chapter  V  were  analyzed  using  principal 
components  analysis,  the  spectral  angle  mapper,  and  a  spectral  matched  filter  to 
determine  whether  thermal  spectral  data  can  be  used  to  detect  targets  in  a  military 
scenario.  The  analysis  process  began  with  an  intensive  study  of  the  Malpai  training 
scene.  This  training  data  was  used  to  establish  a  standard  method  of  applying  the 
hyperspectral  analysis  algorithms. 

A.  MALPAI  COLLECTION  SITE  -  20  OCTOBER  1995 


The  Malpai  site  consisted  of  an  array  of  target  vehicles  and  CCD  materials.  As 
described  in  Chapter  IV,  the  Malpai  targets  were  placed  on  a  surface  of  desert  varnish 
(see  Figure  4.1).  The  background  of  desert  varnish  was  substantially  different  from  the 
targets  placed  on  the  site.  As  a  result,  all  of  the  targets  were  easy  to  detect  throughout  the 
electromagnetic  spectrum.  This  made  the  site  an  excellent  location  for  ground  truth  data 
collection  and  sensor  characterization. 

1.  Data  Processing 

The  Malpai  data  selected  for  analysis  were  collected  at  1000  local  time  on  20 
October.  The  raw  SEBASS  data  were  prepared  using  the  procedures  outlined  in  Chapter 
V,  Section  D.  The  atmospheric  parameters  used  to  convert  the  radiance  measured  by 
SEBASS  to  source  radiance  at  the  ground  were  derived  through  the  use  of  the  plastic 
ruler  algorithm  and  presented  above  in  Figure  5.10  as  examples  of  the  atmospheric 
compensation  process.  During  the  data  preparation  process,  raw  SEBASS  data  was 
converted  into  the  following  data  formats: 
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•  Simulated  FLIR  image  -  single-band  thermal  image 

•  Temperature  image  -  single-band  image  of  estimated  surface  temperature 

•  Radiometrically  calibrated  data  -  radiance  at  sensor 

•  Atmospherically  corrected  data  -  radiance  at  source 

•  Emissivity  data  -  apparent  emissivity 

2.  Single-band  Images 

Two  single-band  images  were  generated  as  outputs  of  the  data  preparation 
process.  The  simulated  FLIR  image  was  created  by  averaging  the  radiance  received  by 
the  sensor  in  each  of  the  128  wavelength  channels.  The  second  single-band  image  was 
composed  of  data  values  representing  the  estimated  surface  temperature  of  each  pixel. 

The  simulated  FLIR  image  of  the  Malpai  scene  is  presented  as  Figure  6.1.  The 
FLIR  image  is  essentially  a  thermal  image  of  the  scene.  The  bright  areas  represent 
materials  emitting  the  most  electromagnetic  energy.  Almost  all  of  the  targets  within  the 
Malpai  collection  site  were  distinguishable  from  the  desert  varnish  because  of  the  high 
temperature  of  the  desert  surface.  Some  man-made  materials,  such  as  the  calibration 
panels  and  several  elements  within  the  JCCD  arrays,  differed  in  temperature  from  the 
other  man-made  objects  and  were  also  easy  to  distinguish.  On  the  other  hand,  the 
vegetation  and  tanks  appeared  to  emit  the  same  amount  of  energy  and  would  be  hard  to 
differentiate  if  they  were  collocated. 

The  single-band  temperature  image  of  the  data  depicted  the  same  information  as 
the  simulated  FLIR  image.  For  comparison,  the  temperature  image  is  presented  as  Figure 
6.2.  Once  again,  almost  all  materials  on  the  training  site  could  be  detected  using 
temperature  differences  alone.  Several  sample  temperatures  have  been  identified  to 
demonstrate  the  thermal  variation  within  the  target  site.  The  dark  desert  varnish  heats  up 
quickly  and  appeared  as  a  bright  surface  in  this  late  morning  image.  It  provided  a  good 
contrast  to  the  cooler  T-72  tanks.  The  high  thermal  inertia  of  the  tanks  caused  them  to 
heat  up  more  slowly.  The  two  US  tankers  in  the  lower  left  comer  were  used  to  determine 
whether  thermal  sensors  can  detect  the  status  of  military  support  vehicles.  The  tanker  on 
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the  left  was  hotter  than  the  tanker  on  the  right  by  8”  K.  The  low  thermal  inertia  of  the 
empty  tanker  caused  it  to  heat  up  quickly  in  the  morning  sun. 

3.  Principal  Components  Analysis 

The  PCA  of  the  Malpai  collection  site  served  two  purposes.  The  analysis  was 
used  to  determine  how  the  results  would  vary  as  the  transform  was  applied  to  the 
radiometrically  calibrated  data,  the  atmospherically  corrected  data,  and  the  emissivity 
data.  It  was  then  used  to  characterize  the  various  targets  within  the  site. 

a.  Radiance  Data  Sets 

The  radiometrically  calibrated  data  and  the  atmospherically  corrected 
data  are  both  in  radiance  units,  sr .  The  calibrated  data  represents  the 

radiance  measured  at  the  sensor  while  the  atmospherically  corrected  data  represents 
radiance  emitted  at  the  surface.  PCA  was  performed  on  both  of  these  radiance  data  sets 
to  determine  if  the  atmospherically  corrected  data  would  provide  substantially  different 
results  from  the  radiometrically  calibrated  data. 

The  results  of  applying  the  principal  components  transformation  to  the  two 
data  sets  are  presented  in  Figures  6.3  and  6.4.  The  figures  portray  how  the  variance  in  the 
original  data  is  represented  in  principal  component  data  space.  Figure  6.3  is  a  graph  of 
the  eigenvalues  of  each  principal  component  (PC)  of  the  radiometrically  calibrated  and 
atmospherically  corrected  data.  The  figure  shows  that  the  transformation  produced 
similar  results  when  applied  to  the  two  data  sets.  The  non-traditional  approach  of  using  a 
logarithmic  horizontal  axis  hides  the  break  in  slope  found  in  log-linear  plots.  In  these 
data,  the  break  occurs  at  PC  band  12.  Figure  6.4  presents  the  same  information  as  a 
percentage  of  original  scene  variance  contained  in  each  PC.  The  figure  shows  that  over 
90%  of  the  variance  within  the  original  image  is  contained  in  the  first  PC.  Each 
subsequent  PC  contains  less  and  less  of  the  information  contained  in  the  original  data. 
After  the  tenth  PC,  the  remaining  components  contain  less  than  0.01  %  of  the  original 
variance  of  the  data. 
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The  first  PC  typically  is  a  weighted  average  of  all  of  the  bands  of  an 
instrument.  In  this  case,  the  first  PC  produces  an  output  similar  to  the  simulated  FT  TP 
image.  Figure  6.5  compares  the  simulated  FLIR  image  to  the  first  PC  of  the  calibrated 
data  and  the  first  PC  of  the  atmospherically  corrected  data.  The  three  images  are  virtually 
identical.  The  fact  that  the  first  PC  contains  90%  of  the  original  variance  in  the  data  and 
represents  the  temperature  variations  within  the  scene  is  a  result  of  the  strong  relationship 
between  temperature  and  emitted  radiation  in  the  LWIR  region  of  the  electromagnetic 
spectrum. 

The  second,  third,  and  fourth  PC  images  of  the  two  types  of  radiance  data 
were  also  similar.  Beginning  with  the  fifth  PC,  the  component  images  of  the  two  data 
sets  began  to  differ.  Figure  6.6  presents  a  comparison  of  the  fifth  PC  image  of  the 
radiometrically  calibrated  data  with  that  of  the  atmospherically  corrected  data.  The  T-72 
tanks  were  easily  detected  in  both  images,  but  were  highlighted  differently  because  of  the 
statistical  differences  between  the  two  data  sets  caused  by  removal  of  atmospheric  effects. 
Despite  the  alteration  in  the  appearance  of  the  individual  images  above  PC  number  four, 
PCA  analysis  of  the  two  data  sets  yielded  the  same  amount  of  information. 

b.  PCA  Results  on  Calibrated  Data  Set 

The  author  chose  to  perform  the  full  PCA  on  the  radiometrically  calibrated 
data  rather  than  the  atmospherically  corrected  data.  Both  data  sets  produced  similar 
results  and  the  use  of  the  calibrated  data  eliminated  any  adverse  effects  a  possible 
mischaracterization  of  the  atmosphere  might  have  on  the  results  of  the  analysis. 

As  mentioned  above,  the  first  PC  of  the  radiometrically  calibrated  data 
was  virtually  identical  in  appearance  to  both  the  simulated  FLIR  image  and  the  surface 
temperature  image.  The  first  PC  image  exploited  the  temperature  variations  within  the 
scene  allowing  materials  to  be  detected  based  on  their  surface  temperature,  thermal  inertia 
properties,  or  internal  heat  sources. 

The  second  PC  provided  a  large  amount  of  information  about  the  scene 
which  was  not  present  in  the  simulated  FLIR  or  first  PC  image.  Figure  6.7  presents  the 
simulated  FLIR  image  alongside  the  second  PC  image  of  the  T-72  array.  The  weighting 
factors  for  this  component  highlighted  materials  which  differed  from  the  desert  varnish 
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background.  As  a  result,  materials  with  high  emissivities  in  the  8.0  to  9.0  )im  region 
where  the  reststrahlen  band  is  located  appeared  bright.  Tank  tracks,  desert  surface 
variations,  and  T-72  array  were  all  clearly  visible. 

The  third  PC  produced  some  additional  discrimination  capabilities  and 
allowed  targets  to  be  differentiated  from  the  background  more  clearly.  A  portion  of  the 
third  PC  image  is  presented  in  Figure  6.8.  The  image  on  the  left  is  the  original  third  PC 
image.  A  brightness  threshold  has  been  applied  to  produce  the  image  on  the  right.  The 
thresholding  highlighted  specific  targets  which  were  emphasized  by  the  transform.  The 
image  clearly  provided  a  means  of  differentiating  between  man-made  and  natural  objects 
within  the  scene.  The  tanks,  decoy  materials,  and  calibration  panels  could  easily  be 
isolated  from  the  background. 

Although  the  second  and  third  PC  images  were  useful  in  discriminating 
between  man-made  objects  and  the  desert  background,  they  did  not  provide  much 
information  to  differentiate  between  different  types  of  man-made  materials.  Several  of 
the  higher  components  provided  this  capability.  Figure  6.9  depicts  the  eighth  PC  image 
with  and  without  thresholding.  The  calibration  panels,  foxhole  cover,  and  bomb  crater 
decoy  were  easily  distinguishable  from  both  the  background  and  the  T-72  tanks. 

Each  of  the  various  PC  images  of  the  radiometrically  calibrated  data 
provided  different  information.  Figure  6.10  displays  how  each  spectral  band  was 
weighted  to  generate  the  first,  second,  and  eighth  PCs.  The  figure  graphically 
demonstrates  that  the  first  PC  is  essentially  a  weighted  average  of  the  spectral  bands.  The 
second  PC  emphasizes  spectral  measurements  within  the  reststrahlen  band  of  the  desert 
sand.  The  determination  of  what  the  weighting  factors  of  the  eighth  PC  represent  is  not 
as  straightforward.  Since  the  eighth  PC  differentiated  between  elements  in  the  JCCD 
array,  the  weighting  function  appears  to  exploit  specific  spectral  differences  between 
these  materials.  Figure  6.1 1  is  a  composite  image  formed  by  combining  these  three  PCs 
into  a  single  red-green-blue  image.  The  figure  characterizes  many  of  the  elements  within 
the  scene.  The  use  of  PC  number  eight  aided  in  the  discrimination  between  the  decoy 
materials,  the  calibration  panels,  and  the  T-72  tanks.  In  the  composite  image,  the  desert 
varnish  is  magenta  and  the  background  vegetation  is  yellow.  Various  man-made  objects 
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vary  in  color  from  orange  in  the  case  of  the  T-72  tanks  to  green  and  blue  in  the  case  of  the 
JCCD  decoy  materials. 

c.  Principal  Components  Analysis  on  Emissivity  Data  Set 

The  principal  components  transformation  was  then  applied  to  the 
emissivity  data  to  determine  how  the  conversion  to  emissivity  would  effect  the  results. 
While  the  eigenvalues  of  the  emissivity  data  followed  the  same  pattern  exhibited  by  the 
eigenvalues  of  the  radiance  data  sets,  the  percentage  of  variance  contained  within  each 
component  data  did  not  drop  below  0.01  %  until  after  the  40*  PC.  Although  this  was  a 
dramatic  difference  from  the  PCA  of  the  radiance  data,  the  PCs  of  the  emissivity  data 
represent  variations  within  the  data  after  the  removal  of  the  temperature  variations  across 
the  scene.  Since  temperature  variations  represent  up  to  90%  of  the  scene  variance, 
removing  this  information  increases  the  proportion  of  remaining  variance  contained  in 
each  PC.  If  the  first  PC  is  removed  from  the  results  of  the  radiometrically  calibrated  or 
atmospherically  corrected  data,  the  variance  graphs  exhibit  a  similar  shape.  This  effect  is 
shown  in  Figure  6. 12.  The  figure  is  a  comparison  of  the  percentage  of  variance  contained 
in  each  PC  of  the  transformed  emissivity  data  and  the  transformed  radiometrically 
calibrated  data  after  the  removal  of  the  first  PC.  The  graph  shows  that  the  first  few  PCs 
of  the  emissivity  data  contain  more  information  than  the  corresponding  PCs  in  the 
calibrated  data.  But,  above  the  tenth  PC,  the  calibrated  data  retains  much  more 
information  within  each  PC. 

Figure  6.13  is  a  shows  the  weighting  functions  used  to  generate  the  first, 
second,  and  fifth  PCs  of  the  emissivity  data.  The  first  component  resembled  the  second 
component  of  the  radiance  data  and  highlighted  materials  which  have  high  emissivities 
within  the  reststrahlen  band  of  the  desert  sand  (see  Figure  6.10).  The  results  of  this 
particular  weighting  are  shown  in  Figure  6.14.  The  higher  components  once  again 
provided  a  means  of  distinguishing  between  various  man-made  materials  within  the 
scene.  The  information  contained  in  the  fourth  and  fifth  PC  images  allowed  for  excellent 
discrimination  between  the  tanks,  the  various  decoy  materials,  and  the  calibration  panels. 
Figure  6.15  is  the  fifth  PC  image.  The  bright  T-72  tanks  were  easy  to  distinguish  from 
the  JCCD  materials  and  calibration  panels. 
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A  color  composite  of  the  second,  fourth  and  fifth  PC  images  is  presented 
as  Figure  6.16.  Within  the  image,  the  T-72s  appear  yellow,  the  vegetation  is  orange,  and 
elements  of  the  JCCD  array  are  blue-green.  Tank  tracks  are  highlighted  in  dark  blue  or 
magenta.  Because  this  image  is  based  entirely  on  emissivity  data,  it  has  no  thermal 
component  and  would  produce  the  same  results  regardless  of  illumination  conditions. 
This  image  characterized  the  scene  well,  but  the  composite  image  generated  from  the 
radiometrically  calibrated  data  (Figure  6. 1 1)  contained  more  information.  The  tanks 
tracks,  vegetation,  and  JCCD  materials  were  more  distinct  in  Figure  6.1 1.  The 
temperature  component  of  the  radiance  data  provided  some  additional  discrimination 
capabilities. 

4.  Spectral  Angle  Mapper  Analysis 

The  spectral  angle  mapper  was  applied  to  the  emissivity  data  generated  from  the 
Malpai  collection  site.  As  described  in  Chapter  E,  Section  C,  the  SAM  compares  the 
emissive  spectra  within  a  data  set  to  a  library  of  reference  spectra.  The  Malpai  training 
set  provided  an  initial  opportunity  to  see  how  well  the  ground  truth  spectra  matched  the 
data  collected  by  SEBASS.  It  also  served  as  the  starting  point  for  the  development  of  a 
spectral  library  which  would  be  used  on  the  data  collected  over  the  JCCD  and  Arroyo 
sites.  The  ground  truth  spectra  collected  by  John  Hopkins  University  served  as  the 
primary  source  of  reference  spectra.  The  library  was  augmented  with  spectra  extracted 
from  the  Malpai  Collection  Site  data. 

a.  SAM  Analysis  Using  Ground  Truth  Spectra 

The  initial  SAM  analysis  of  the  Malpai  collection  site  utilized  only  ground 
truth  spectra  collected  by  John  Hopkins  University  (see  Table  5.1).  The  ground  truth  data 
included  two  target  vehicle  spectra.  The  T1326_JHU  spectrum  was  collected  from  the 
center  T-72  in  the  five  tank  array.  The  MBTl  152_JHU  spectrum  was  collected  from  one 
of  the  T-80  surrogate  tanks. 

The  first  SAM  trial  demonstrated  that  the  ground  truth  spectra  provided  a 
good  match  to  some  of  the  materials  in  the  scene.  The  algorithm  was  run  with  a  spectral 
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angle  threshold  of  0.03  radians.  The  results  are  presented  as  Figure  6. 17.  All  of  the 
military  vehicles  were  detected  as  one  of  the  two  target  vehicle  types  and  displayed  as 
either  red  or  blue.  The  algorithm  also  detected  many  of  the  tracks  left  by  the  vehicles  as 
they  were  moved  into  position.  The  algorithm  did  not  correctly  identify  the  vegetation, 
some  of  the  lighter  sand,  and  portions  of  the  tank  tracks.  The  vegetation  was 
misclassified  as  Air  Force  desert  camouflage  and  displayed  as  orange.  Portions  of  the 
vehicle  tracks  were  misclassified  as  T-80  targets  and  portrayed  as  blue.  Various  JCCD 
materials  were  misclassified  due  to  the  lack  of  appropriate  reference  spectra.  These 
materials  were  prototype  materials  and  did  not  correspond  well  to  the  ground  truth 
camouflage  materials. 

b.  SAM  Analysis  Using  Ground  Truth  and  Extracted  Spectra 

To  improve  classification  accuracy,  the  reference  library  was  expanded  to 
include  spectra  extracted  from  the  Malpai  data.  Spectra  representing  the  desert  varnish, 
the  lighter  desert  sand,  vegetation,  and  a  disturbed  desert  surface  resulting  from  vehicular 
traffic  were  extracted  from  known  locations  within  the  Malpai  scene.  Since,  the  foxhole 
camouflage  and  IR  applique  exhibited  very  distinctive  spectral  characteristics  and  did  not 
correspond  to  any  of  the  ground  truth  camouflage  materials,  they  were  also  added  to  the 
reference  library.  The  full  reference  library  consisted  of  the  materials  listed  in  Tables  5.1 
and5.2. 

The  use  of  the  expanded  reference  library  improved  the  results  of  the  SAM 
analysis  dramatically.  The  vegetation,  variations  in  desert  background,  and  JCCD 
materials  were  classified  with  much  greater  accuracy.  The  false  alarms  involving 
misclassification  of  the  desert  background  and  vegetation  were  also  eliminated.  Figure 
6.18  is  a  portion  of  the  resulting  classification  image  using  a  spectral  angle  threshold  of 
0.03  radians.  The  tank  tracks,  portrayed  in  orange,  were  clearly  distinct  from  the  lighter 
sand  and  desert  varnish.  The  vegetation  was  also  clearly  identified.  But,  the 
improvement  in  background  classification  and  the  reduction  in  false  alarms  was  achieved 
at  the  cost  of  easily  locating  all  of  the  threat  vehicles.  Only  four  of  the  T-72s  were 
positively  classified  as  being  either  T-72  or  T-80  targets.  The  remaining  T-72  had  only 
seven  pixels  correctly  classified  as  matching  the  T-72  spectrum. 
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5.  Spectral  Matched  Filter  Analysis 


The  SMF  identifies  spectra  within  a  scene  that  closely  match  an  input  reference 
spectrum.  The  algorithm  was  applied  to  the  Malpai  data  to  determine  how  well  each  of 
the  target  T-72s  matched  the  ground  tmth  spectra.  It  was  also  used  to  conduct  a  quick 
terrain  classification  of  the  scene.  Reference  spectra  for  the  SMF  were  taken  from  the 
same  reference  library  generated  for  the  SAM  analysis. 

The  SMF  was  first  employed  to  detect  the  military  target  vehicles  using  the  T-72 
and  T-80  spectra.  The  slight  difference  between  the  T-72  and  T-80  spectra  was  evident 
in  the  SMF  results.  Figure  6. 19  is  a  portion  of  the  SMF  results  for  the  Malpai  data.  The 
colors  have  been  reversed  in  the  figure  to  better  highlight  the  SMF  results.  The  dark 
areas  represent  spectra  which  match  the  T-72  and  T-80  reference  spectra.  The  image  on 
the  left  is  the  product  of  a  matched  filter  using  the  T-72  spectra.  Only  three  of  the  five 
vehicles  are  positively  detected,  but  the  vehicles  that  were  detected  were  well  defined  and 
distinct  from  the  background.  The  image  on  the  right  is  the  result  of  the  T-80  trail.  All  of 
the  T-72  tanks  were  detected  within  the  image,  but  the  signal-to-clutter  ratio  was  much 
lower  making  them  harder  to  distinguish  from  the  background.  Although  not  portrayed  in 
the  figure,  the  two  T-80  targets  at  the  bottom  of  the  collection  site  were  strongly 
highlighted  in  T-80  filtered  image.  In  both  images,  the  vegetation  was  not  highlighted 
indicating  that  the  spectra  from  these  targets  were  sufficiently  distinct  to  distinguish  them 
from  vegetation.  When  reference  spectra  were  extracted  from  the  CCD  materials  at  the 
north  end  of  the  site,  the  SMF  accurately  classified  the  materials. 

Terrain  classification  using  the  SMF  also  worked  well.  The  vegetation  within  the 
washes  was  accurately  classified  using  the  extracted  vegetation  spectra.  The  use  of  a 
single  desert  varnish  spectrum  allowed  several  desert  features  to  be  characterized 
according  to  how  much  they  deviated  from  the  reference  spectmm.  In  fact,  the  use  of  the 
desert  varnish  spectmm  produced  a  more  extensive  vehicle  traffic  assessment  than  use  of 
the  spectmm  extracted  directly  from  a  particular  tank  track.  Figure  6.20  displays  the 
results  of  application  of  the  SMF  with  desert  varnish  as  the  reference  spectmm.  The  tank 
tracks  are  dark  because  they  differ  from  the  desert  varnish  which  has  a  deeper  reststrahlen 
emissivity  minimum. 
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B.  JCCD  COLLECTION  SITE  -  23  OCTOBER  1995 


The  JCCD  collection  site  consisted  of  three  rows  of  vehicles  as  described  in 
Chapter  IV,  Section  B.  Two  of  the  rows  contained  decoy  targets  and  a  third  row 
contained  actual  threat  vehicles.  The  SEBASS  collection  mission  conducted  at  1000 
local  time  on  23  October  included  flight  lines  directly  over  all  three  rows. 

1.  Data  Processing 

The  SEBASS  data  was  prepared  using  the  same  procedures  outlined  in  Chapter  V, 
Section  D.  Atmospheric  parameters  were  calculated  for  each  of  the  three  passes.  The 
collection  runs  over  the  Vehicle  Target  Line  contained  much  more  vegetation  than  the 
passes  over  the  two  decoy  lines.  As  a  result,  the  atmospheric  parameters  derived  from  the 
Vehicle  Target  Line  were  much  more  reliable  and  compared  more  favorably  with  the 
ARTSPEC  model.  Since  all  three  data  sets  were  collected  in  a  span  of  20  minutes  over 
the  same  area,  the  parameters  derived  from  the  Vehicle  Target  Line  data  were  applied  to 
all  three  data  sets.  The  atmospheric  parameters  are  presented  in  Figure  6.21.  The 
ARTSPEC  model  results  are  depicted  as  a  dotted  line  and  the  scaled  plastic  ruler  results 
are  presented  as  a  solid  line. 

2.  JCCD  Decoy  Line  One  -  23  October  1995 

The  first  JCCD  decoy  line  consisted  of  five  decoys.  The  targets  included  murals 
depicting  T-72  and  M-1  tanks  as  well  as  decoys  constructed  of  scrap  materials.  The  CCD 
preparations  of  the  five  decoys  were  described  in  Figure  4.2  and  Table  4.2. 

a.  Single-band  Images 

The  simulated  FLIR  image  of  Decoy  Line  One  is  presented  as  Figure  6.22. 
Decoy  One  is  easily  detectable  in  the  FLIR  image  due  to  its  rectangular  structure  and  the 
hot  spot  directly  in  the  center  of  the  target.  Its  shape  and  thermal  properties  also  allow  the 
decoy  to  be  differentiated  from  a  true  tank.  Figure  6.23  is  a  comparison  of  the  FLIR 
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image  of  Decoy  One  with  a  FLffi.  image  of  a  T-72  from  the  Malpai  data.  Spatially  the 
two  targets  appear  very  distinct.  In  addition,  the  temperature  of  the  decoy  varied  from 
308°  K  in  the  center  to  299°  K  on  the  outer  edges.  Although  the  real  T-72s  also  exhibited 
a  temperature  variation  between  the  turret  and  body,  the  turrets  were  cooler  than  the  body 
during  this  time  of  day  due  to  the  thermal  inertia  of  their  armor  plating 

b.  Principal  Components  Analysis 

Principal  components  analysis  was  performed  on  the  radiometrically 
calibrated  data  and  apparent  emissivity  data.  All  five  decoys  were  detected  during  the 
PCA  analysis  of  the  radiometrically  calibrated  data.  Analysis  of  the  emissivity  data  also 
resulted  in  the  detection  of  the  five  decoys. 

The  results  of  the  PCA  of  the  radiometrically  calibrated  were  similar  to  the 
results  of  PCA  of  the  Malpai  data.  The  first  PC  image  was  identical  in  appearance  to  the 
FLIR  and  surface  temperature  images  and  represented  temperature  variations  within  the 
scene.  The  second  PC  highlighted  tank  tracks,  vegetation,  and  targets  which  differed 
from  the  majority  background  which  in  this  case  was  the  light  sand.  This  PC  image  is 
presented  as  Figure  6.24.  Decoy  One,  Decoy  Four,  and  Decoy  Five  appeared  as  bright 
objects.  In  addition,  various  tracks  and  trails  provided  an  indication  of  vehicular  traffic 
and  served  as  clues  to  the  possible  location  of  other  targets  and  decoys  within  the  scene. 
Although  Decoy  Three  was  not  directly  identifiable,  the  edge  of  a  camouflage  net  was 
clearly  evident.  Thresholding  was  applied  to  an  enlarged  portion  of  the  image  to  more 
clearly  differentiate  the  edge  of  the  camouflage  netting  from  the  general  background.  The 
higher  order  PC  images  provided  additional  information  about  the  decoys. 

A  color  composite  of  the  first  three  PC  images  was  used  to  combine  the 
information  contained  in  three  components.  The  first  three  PC  images  were  selected  to 
characterize  the  desert  background  and  highlight  any  man-made  objects.  The  image  is 
presented  as  Figure  6.25.  Four  of  the  decoys  could  be  detected  and  the  general  JCCD 
background  was  characterized.  Although  Decoy  Two  was  detected  in  the  third  PC  image, 
its  small  size  made  it  difficult  to  detect  in  the  composite  image.  In  the  figure,  the  desert 
varnish  appears  as  magenta.  Light  sand,  disturbed  sand,  and  tank  tracks  appear  as  various 
shades  of  blue-green.  The  vegetation  is  light  green-orange.  Decoy  One  is  bright  green. 
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The  other  four  decoys  are  yellow- white.  Table  6.1  summarizes  the  color  index  produced 
by  this  particular  combination  of  PC  images. 


Color 

Material 

Yellow-White 

Decoy/Target 

Green-Orange 

Vegetation 

Magenta-Red 

Desert  Surface  (Desert  Varnish) 

Blue 

Tank  Tracks 

Table  6.1.  Color  Index  for  Principal  Component  Image  Using  PC  Bands  One,  Two,  and 
Three  from  Radiometrically  Calibrated  Data. 


c.  Spectral  Angle  Mapper  Analysis 

The  spectral  angle  mapper  was  applied  to  the  emissivity  using  the 
reference  library  generated  during  the  SAM  analysis  of  the  Malpai  data.  The  first  run 
employed  a  spectral  threshold  of  0.03  radians.  The  resulting  classification  image  is 
presented  as  Figure  6.26.  The  algorithm  detected  four  of  the  five  decoys.  Decoy  One 
highlighted  in  black  because  it  was  constructed  out  of  aluminum  matting  which  has  a 
spectrum  which  does  not  resemble  any  of  the  input  spectra.  The  cyan  cluster  above 
Decoy  One  is  a  ghost  image  of  the  decoy  caused  by  the  reflection  and  diffraction  artifact 
identified  in  Chapter  V,  Section  B.  Decoy  Three  and  Decoy  Five  employed  various  types 
of  woodland  camouflage  and  were  classified  as  LCSS  woodland  camouflage.  Decoy 
Four  was  classified  as  similar  to  the  prototype  foxhole  camouflage.  The  decoy  was  a 
visual  mural  of  a  T-72  placed  on  top  of  woodland  camouflage.  No  matches  were  found 
for  T-72  or  T-80  spectra.  The  algorithm  had  some  problems  classifying  the  background. 
It  misclassified  some  of  the  tank  trails  as  vegetation.  Decreasing  the  spectral  angle 
threshold  to  0.01  radians  produced  similar  results. 

To  determine  the  extent  of  the  misclassification  problem  involving  the 
vegetation  reference  spectmm  and  the  tank  tracks,  a  SAM  trial  was  run  with  the  extracted 
vegetation  spectrum  excluded.  This  second  trial  correctly  identified  the  vehicle  tracks 
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and  light  sand.  It  revealed  that  much  of  the  sand  within  the  wash  had  been  classified  as 
vegetation. 

d.  Spectral  Matched  Filter  Analysis 

The  SMF  analysis  of  Decoy  Line  One  focused  on  detecting  the  various 
targets  within  the  scene.  The  algorithm  was  initially  run  with  the  1-12  and  T-80  spectra 
as  inputs  to  determine  if  any  targets  could  be  quickly  classified  as  threat  vehicles.  The 
various  camouflage  schemes  were  then  used  as  reference  spectra  for  detection  of  any 
CCD  efforts.  Finally,  the  desert  varnish  spectra  was  used  to  determine  whether  any 
vehicle  tracks  could  be  identified. 

The  T-72  and  T-80  spectra  reference  produced  detections  of  Decoy  One 
and  Decoy  Four.  Both  of  these  detections  had  low  signal-to-clutter  ratios  and  were  based 
in  large  part  on  the  large  difference  these  spectra  had  from  the  background.  The  LOSS 
woodland  camouflage  spectrum  produced  detections  of  all  of  the  decoys  except  Decoy 
Five.  The  resulting  filtered  image  is  presented  as  Figure  6.27.  To  increase  the  contrast  of 
the  detected  targets,  high  response  values  have  been  coded  as  dark  gray.  Decoy  One  was 
highlighted  even  though  it  was  not  camouflaged.  Decoy  Two  was  covered  with  a  LOSS 
desert  camouflage  and  incorrectly  classified  as  a  woodland  camouflage  match.  Decoy 
Three  employed  an  AFCENT  woodland  camouflage  which  appeared  to  be  similar  to  the 
Less  woodland  camouflage.  Decoy  Four  employed  the  LOSS  woodland  camouflage  and 
had  very  high  signal-to-clutter  ratio.  The  use  of  the  German  woodland  camouflage 
spectrum  as  a  reference  spectrum  produced  similar  results.  When  the  SMF  was  applied 
with  the  prototype  foxhole  camouflage  as  the  reference,  the  algorithm  detected  Decoy 
Four  with  a  very  high  signal-to-clutter  ratio  indicating  that  the  decoy  was  a  good  match 
for  the  spectrum. 

The  application  of  the  tank  track  spectrum  and  the  vegetation  spectrum 
produced  good  ground  characterization  results.  When  the  extracted  vegetation  spectrum 
was  used  in  the  SMF  algorithm,  all  of  the  decoys  near  vegetation  were  easily  detected 
because  they  were  marked  by  dark  patches  within  the  bright  vegetation.  Use  of  the  tank 
track  spectrum  revealed  the  tracks  leading  to  Decoy  Five.  This  was  the  only  indication  of 
the  fifth  decoy  produced  by  the  SMF. 
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3.  JCCD  Decoy  Line  Two  -  23  October  1995 


The  second  JCCD  decoy  line  consisted  of  four  decoy  targets.  The  targets  included 
two  high  fidelity  mobile  artillery  mockups,  a  field  expedient  decoy  constructed  of  scrap 
material,  and  a  mural  of  an  M-1  tank.  Figure  4.2  and  Table  4.2  describe  the  CCD 
preparations  of  the  four  decoys. 

a.  Single-band  Images 

The  simulated  FLIR  image  of  JCCD  Decoy  Line  Two  is  presented  as 
Figure  6.28.  Decoy  One  is  easily  detectable  in  the  FLIR  image  due  to  its  rectangular 
structure  and  the  hot  spot  within  the  camouflage.  The  target’s  linear  shape  differentiates 
it  from  the  rounded,  irregular  vegetation  in  the  wash.  As  was  the  case  with  the  previous 
data,  the  surface  temperature  image  (not  shown)  was  similar  to  the  simulated  FLIR 
image. 

b.  Principal  Components  Analysis 

Principal  components  analysis  was  performed  on  both  the  radiometrically 
calibrated  data  and  apparent  emissivity  data.  The  PCA  analysis  of  the  radiometrically 
calibrated  data  resulted  in  the  detection  of  Decoys  One,  Decoy  Two,  and  Decoy  Four. 
Analysis  of  the  emissivity  data  resulted  in  the  detection  of  the  all  four  decoys. 

The  PCA  results  of  the  radiometrically  calibrated  data  were  similar  to  the 
results  of  the  Malpai  data.  The  second  PC  highlighted  tank  tracks,  vegetation,  and  targets 
which  differed  from  the  desert  sand.  Figure  6.29  demonstrates  how  the  second  PC  image 
highlighted  Decoy  One  and  Decoy  Four.  The  portion  of  the  image  in  which  Decoy  Four 
was  located  has  been  enlarged  to  further  demonstrate  the  detection  capabilities  of  this  PC. 
The  third  PC  image  highlighted  Decoy  Two  and  Decoy  Four. 

A  color  composite  of  the  first  three  PC  images  is  presented  as  Figure  6.30. 
This  image  illustrates  the  detection  of  three  of  the  decoys.  The  resulting  image  is  similar 
to  the  three-band  composite  of  Decoy  Line  One  presented  as  Figure  6.25.  The  color 
index  is  similar  to  that  listed  in  Table  6.1.  The  desert  background  is  magenta.  Tracks  and 
disturbed  sand  are  blue.  The  vegetation  is  a  light  orange-green.  Decoy  One  is  also  green. 
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but  has  a  darker  hue  and  a  rectangular  outline.  Decoy  Two  and  Decoy  Four  are  yellow- 
white.  Decoy  Three,  although  not  directly  detectable,  may  be  located  at  the  base  of  the 
tracks  highlighted  in  the  image  on  the  right. 

The  analysis  of  the  emissivity  data  provided  an  opportunity  to  detect 
Decoy  Three.  The  decoy  was  detected  in  the  third  and  fourth  PCs.  Figure  6.31  is  a  subset 
of  a  color  composite  of  the  first  three  PC  images  of  the  emissivity  data.  It  is  a  portion  of 
the  decoy  line  containing  Decoy  Three  and  Decoy  Four.  Although  the  background  colors 
within  this  image  differ  from  those  in  Figure  6.30,  the  decoys  still  appear  as  white. 

c.  Spectral  Angle  Mapper  Analysis 

The  spectral  angle  mapper  was  applied  to  the  emissivity  data  using  the 
same  reference  library  generated  during  the  SAM  analysis  of  the  Malpai  data.  The  first 
run  employed  a  spectral  angle  threshold  of  0.03  radians.  Only  Decoy  Four  was  identified. 
Decoy  Four  was  a  mural  of  an  M-1  with  an  East  German  woodland  camouflage  net  in 
front.  The  algorithm  classified  it  as  a  combination  of  LCSS  woodland  camouflage  and 
Air  Force  desert  camouflage.  The  slight  misclassification  was  due  to  the  lack  of  ground 
truth  on  the  East  German  camouflage.  Once  again,  much  of  the  desert  sand  including 
some  of  the  tank  tracks  was  incorrectly  classified  as  vegetation.  The  algorithm  produced 
similar  results  when  run  with  a  smaller  spectral  angle  threshold  of  0.01  radians. 

The  analysis  was  repeated  with  a  spectral  library  that  excluded  the 
extracted  vegetation  spectrum.  The  results  are  presented  in  Figure  6.32.  In  order  to 
separate  the  T-72/T-80  targets  from  the  Air  Force  desert  camouflage,  the  Air  Force  desert 
camouflage  is  portrayed  as  black  in  this  image.  Despite  a  few  false  alarms.  Decoy  One, 
Decoy  Two,  and  Decoy  Four  were  detected.  Decoy  One  was  correctly  identified  as  being 
covered  with  Air  Force  desert  camouflage.  Decoy  Two  was  classified  as  a  combination 
of  T-72  and  Air  Force  desert  camouflage  spectra.  It  was  actually  a  mobile  artillery  decoy 
covered  by  a  woodland  camouflage  net.  Decoy  Three  was  not  clearly  detected  even 
though  the  German  camouflage  employed  on  the  target  was  one  of  the  ground  truth 
spectra  used  in  the  reference  library.  Decoy  Four  was  classified  as  a  combination  of  T- 
72,  Air  Force  desert  camouflage,  and  LCSS  woodland  camouflage.  It  was  actually  a 
mural  of  an  M-1  tank  next  to  an  East  German  camouflage  net.  A  spatial  filter  was  then 
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applied  to  remove  some  of  the  false  alarms.  Although  this  reduced  the  number  of  small 
false  alarms,  several  tank-sized  false  alarms  still  remained.  An  example  of  the  spatially 
filtered  output  displaying  Decoy  Two,  Decoy  Four,  and  several  false  alarms  is  presented 
as  Figure  6.33. 

d.  Spectral  Matched  Filter  Analysis 

The  SMF  analysis  of  JCCD  Decoy  Line  Two  detected  all  four  decoys 
using  several  of  the  ground  truth  camouflage  spectra.  Decoy  Two,  Decoy  Three,  and 
Decoy  Four  were  detected  using  the  LCSS  and  West  German  woodland  camouflage 
spectra.  The  resulting  filtered  image  is  presented  as  Figure  6.34.  Dark  pixels  represent 
data  which  closely  match  the  input  reference  spectrum.  All  three  of  these  decoys 
employed  some  form  of  woodland  camouflage.  In  a  separate  application  of  the  SMF  with 
the  Air  Force  desert  camouflage  spectra,  the  large  net  over  Decoy  One  was  classified. 
Application  of  the  T-72  and  T-80  spectra  did  not  produce  any  significant  target 
detections. 


4.  JCCD  Vehicle  Target  Line  -  23  October  1995 

The  final  JCCD  line  consisted  of  the  five  threat  vehicles.  The  targets  included 
four  T-72  tanks  and  one  mobile  artillery  piece.  Figure  4.2  and  Table  4.3  describe  the 
CCD  preparations  of  the  targets.  These  targets  were  camouflaged  using  earth,  sand  bags, 
and  cut  vegetation.  Only  one  of  the  targets  employed  camouflage  material  as  its  primary 
CCD  effort. 

a.  Single-band  Images 

The  simulated  FLIR  image  of  JCCD  Vehicle  Target  Line  is  presented  in 
Figure  6.35.  The  target  line  was  a  wash  with  a  moderate  amount  of  vegetation.  The 
SEE  ASS  target  run  included  Tanks  One  through  Four.  The  engine  of  Tank  Four  was 
operated  during  the  morning.  The  heat  generated  in  the  engine  allowed  the  FLIR  and 
surface  temperature  images  to  detect  the  tank.  The  temperature  of  the  hottest  part  of  the 
tank  was  325°  K.  No  other  vehicles  were  detected  in  the  single-band  images. 
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b.  Principal  Components  Analysis 


The  PCA  performed  on  the  radiometrically  calibrated  data  and  apparent 
emissivity  data  produced  similar  results.  Only  one  target  was  detected.  Tank  One  was 
first  detected  in  the  third  PC  image  of  the  radiometrically  calibrated  data.  It  was  detected 
in  the  second  PC  of  the  emissivity  data. 

A  color  composite  of  the  first  three  PC  images  is  presented  as  Figure  6.36. 
This  image  highlights  Tank  One  and  characterizes  of  the  general  JCCD  background.  The 
colors  for  this  image  are  similar  to  those  listed  in  Table  6.1.  The  desert  background  is 
magenta.  Tracks  and  disturbed  sand  are  blue.  The  vegetation  in  this  image  is  a  more 
distinct  green.  The  one  detected  target  is  yellow-white. 

c.  Spectral  Angle  Mapper  Analysis 

The  spectral  angle  mapper  was  applied  to  the  emissivity  data  using  the 
reference  library  generated  during  the  SAM  analysis  of  the  Malpai  data.  The  first  run 
employed  a  spectral  angle  threshold  of  0.03  radians.  The  algorithm  detected  the  desert 
camouflage  over  Tank  One  as  being  similar  to  Air  Force  desert  camouflage.  The  vehicle 
was  actually  camouflaged  with  a  United  Kingdom  desert  camouflage  pattern.  The 
algorithm  also  detected  a  set  of  tank  tracks  that  passed  through  the  wash.  Figure  6.37  is 
the  classification  image  generated  after  the  algorithm  was  run  with  a  smaller  spectral 
angle  threshold. 

d.  Spectral  Matched  Filter  Analysis 

The  spectral  matched  filter  analysis  of  the  JCCD  Target  Vehicle  Line 
produced  one  detection  based  on  CCD  efforts.  No  T-72  or  T-80  targets  were  detected. 

The  camouflage  employed  on  Tank  One  was  identified  in  several  filtered  images.  The 
two  forms  of  desert  camouflage.  Air  Force  desert  camouflage  and  LCSS  desert 
camouflage,  detected  the  United  Kingdom  desert  camouflage  over  the  target.  The 
prototype  foxhole  camouflage  spectrum  also  produced  a  good  detection.  Of  the  three 
camouflage  spectra  that  detected  the  target,  the  Air  Force  desert  camouflage  filter 
achieved  the  highest  signal-to-clutter  ratio  indicating  that  the  camouflage  employed  by 


67 


the  tank  resembled  it  the  closest.  The  results  of  the  SMF  analysis  using  the  Air  Force 
camouflage  are  presented  in  Figure  6.38.  The  image  is  scaled  such  that  darker  pixels 
represent  matches  to  the  target  spectram.  The  figure  illustrates  the  detection  of  Tank 
One. 


C.  ARROYO  COLLECTION  SITE  -  23  OCTOBER  1995 

The  Arroyo  site  contained  two  TBM  targets  and  four  threat  vehicles.  The  threat 
vehicles  were  arranged  in  a  protective  perimeter  around  the  high  value  targets.  The  site 
layout  is  depicted  in  Figure  4.3.  The  CCD  measures  employed  on  each  target  are  listed  in 
Table  4.4. 


1.  Data  Processing 

The  data  collected  over  the  Arroyo  Collection  Site  was  prepared  using  the  same 
procedures  outlined  in  Chapter  V,  Section  D.  The  Arroyo  site  had  a  moderate  amount  of 
vegetation.  As  a  result,  the  plastic  ruler  atmospheric  compensation  method  performed 
well.  The  scaled  atmospheric  parameters  and  results  of  the  ARTSPEC  model  are 
presented  in  Figure  6.39. 

2.  Single-band  Images 

The  simulated  FLIR  image  of  Arroyo  collection  site  is  presented  as  Figure  6.40. 
The  SEBASS  collection  mission  flew  directly  over  the  two  TBM  assets.  One  of  the 
TBMs  was  visible  in  the  both  the  FLIR  image  and  the  surface  temperature  image.  The 
large  shape  of  TBM  One  in  relation  to  the  surrounding  vegetation  made  it  a  clear  target. 
A  structure  just  north  of  the  road  in  the  bottom  of  the  image  is  another  possible  target.  It 
corresponds  to  the  location  of  Tank  Four  which  was  covered  with  a  thermal  space  frame. 
The  other  targets  in  the  scene  were  not  be  detected  in  either  the  FT  JR  or  surface 
temperature  image. 
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3.  Principal  Components  Analysis 


Principal  components  analysis  of  the  radiometrically  calibrated  data  and 
emissivity  data  resulted  in  the  detection  of  the  two  TBMs  and  Tank  Four.  The  first 
component  of  the  radiance  data  resembled  the  FLIR  and  temperature  images.  The  second 
highlighted  vegetation,  tracks  within  the  sand,  and  man-made  objects.  The  PC  images 
provided  some  interesting  information  concerning  TBM  One.  The  camouflage  over  TBM 
One  was  composed  of  desert  camouflage  with  some  woodland  camouflage  patches. 

These  patches  were  visible  in  the  different  principal  component  images.  Figure  6.41 
shows  a  portion  of  the  second  and  third  principal  component  images  in  which  TBM  One 
is  visible.  The  variations  within  the  camouflage  netting  as  well  as  the  boundary  between 
the  camouflage  and  the  desert  varnish  are  clearly  visible.  Analysis  of  the  emissivity  data 
resulted  in  detection  of  the  same  three  targets. 

4.  Spectral  Angle  Mapper  Analysis 

The  spectral  angle  mapper  was  applied  to  the  emissivity  data  using  the  same 
reference  library  generated  during  the  SAM  analysis  of  the  Malpai  data.  The  resulting 
classification  image  of  this  trial  is  presented  as  Figure  6.42.  A  spectral  angle  threshold  of 
0.03  radians  was  employed.  The  algorithm  identified  the  two  TBMs  and  Tank  Four. 

Tank  Four  and  TBM  Two  were  classified  as  resembling  the  prototype  foxhole 
camouflage  and  are  portrayed  in  purple.  Examination  of  the  spectra  of  these  targets 
revealed  that  their  emissive  spectra  matched  that  of  the  foxhole  camouflage  very  closely. 
TBM  One  was  classified  as  employing  a  combination  of  desert  and  woodland 
camouflage.  Various  tracks  and  trails  were  identified  and  indicated  that  a  lot  of  vehicle 
traffic  had  occurred  on  the  site.  Several  LOSS  woodland  camouflage  false  alarms  were 
noted  within  some  of  the  vegetation.  In  addition,  an  area  of  disturbed  desert  sand  in  the 
lower  right  comer  of  the  figure  was  misclassified.  The  algorithm  classified  it  as  a  T- 
72/T-80  target.  Examination  of  the  spectra  in  this  area  revealed  that  it  did  not  match  the 
other  desert  surface  spectra  very  well.  Therefore,  the  SAM  algorithm  chose  the  closest 
spectra  which  was  the  T-72  target  spectra. 
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Based  on  the  results  of  the  JCCD  analysis,  the  SAM  was  applied  once  again 
without  the  extracted  vegetation  spectrum  to  see  the  results  of  any  possible 
misclassification  of  the  desert  sand.  This  second  trial  reduced  the  amount  of  desert  sand 
within  die  wash  that  was  classified  as  vegetation.  It  also  revealed  that  some  of  the 
woodland  camouflage  over  TBM  One  had  been  classified  as  vegetation.  The  same  three 
targets  were  detected,  but  TBM  One  identified  with  a  higher  confidence  level. 

5.  Spectral  Matched  Filter  Analysis 

Application  of  the  SMF  resulted  in  the  characterization  of  both  TBM 
emplacements  and  Tank  Four.  The  combination  woodland  and  desert  camouflage  of 
TBM  One  was  detected  using  the  LOSS  woodland  and  desert  camouflages.  Figure  6.43  is 
the  resulting  filtered  image.  The  darker  pixels  represent  matches  to  the  target  spectra. 

The  difference  in  appearance  of  the  camouflage  between  the  two  filtered  images  was 
caused  by  use  of  patches  of  woodland  camouflage  within  the  desert  camouflage.  Tank 
Four  and  TBM  Two  were  easily  detected  using  the  prototype  foxhole  camouflage.  The 
camouflage  scheme  employed  by  these  two  targets  appears  to  be  based  on  the  same 
chemical  composition. 

D.  TARGET  DETECTION  RESULTS 

The  results  of  the  various  algorithms  are  presented  below  in  Table  6.2.  The 
targets  and  decoys  are  labeled  as  detected  if  they  were  picked  out  of  the  background  by 
the  algorithm.  A  target  or  decoy  is  labeled  as  characterized  if  specific  materials  or 
characteristics  about  the  target  were  identified.  This  included  cases  where  a  material  so 
closely  matched  a  particular  ground  truth  spectra  that  the  materials  must  be  related.  For 
example,  no  ground  truth  was  available  for  the  camouflage  materials  employed  over 
TBM  Two  and  Tank  Four  in  the  Arroyo  site.  The  spectra  of  these  materials  matched  that 
of  the  prototype  foxhole  camouflage  very  closely  and  indicated  that  the  materials  must  be 
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related.  In  this  case,  a  match  between  the  foxhole  camouflage  and  the  CCD  material  was 
considered  a  characterization. 


Target 

FUR 

PCA 

(radiometrically 
calibrated  data) 

PCA  (emissivity 
Data) 

SAM 

SMF 

JCCD  Line  1 

Decoy  1 

Detected 

Detected 

Detected 

Detected 

Detected 

Decoy  2 

Not  Detected 

Detected 

Detected 

Not  Detected 

Detected 

Decoy  3 

Not  Detected 

Detected 

Detected 

Characterized 

Characterized 

Decoy  4 

Not  Detected 

Detected 

Detected 

Characterized 

Characterized 

Decoys 

Not  Detected 

Detected 

Detected 

Characterized 

Not  Detected 

JCCD  Line  2 

Decoy  1 

Detected 

Detected 

Detected 

Characterized 

Characterized 

Decoy  2 

Not  Detected 

Detected 

Detected 

Detected 

Characterized 

Decoy  3 

Not  Detected 

Not  Detected 

Detected 

Not  Detected 

Characterized 

Decoy  4 

Not  Detected 

Detected 

Detected 

Characterized 

Characterized 

JCCD  Target  Vehicle  Line 

Tank  1 

Not  Detected 

Detected 

Detected 

Characterized 

Characterized 

Tank  2 

Not  Detected 

Not  Detected 

Not  Detected 

Not  Detected 

Not  Detected 

Tank  3 

Not  Detected 

Not  Detected 

Not  Detected 

Not  Detected 

Not  Detected 

Tank  4 

Detected 

Not  Detected 

Not  Detected 

Not  Detected 

Not  Detected 

Arroyo 

TBM1 

Detected 

Detected 

Detected 

Characterized 

Characterized 

TBM2 

Not  Detected 

Detected 

Detected 

Characterized 

Characterized 

Tank  4 

Detected 

Detected 

Detected 

Characterized 

Characterized 

Table  6.2.  Results  of  Hyperspectral  Analysis  of  SEE  ASS  Data  Collected  20  October  1995 


at  1000  and  23  October  1995  at  1000. 
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Vn.  SUMMARY  AND  CONCLUSIONS 


The  first  operations  of  a  thermal  infrared  imaging  spectrometer  were  conducted  on 
20  -  24  October  1995.  Analysis  of  the  LWIR  spectral  data  collected  resulted  in  the 
successful  detection  of  military  targets.  Targets  were  detected  through  principal 
components  analysis  of  the  radiometrically  calibrated  SEBASS  data.  Analysis  of 
processed  SEBASS  resulted  in  the  classification  and  characterization  of  targets  and 
background  materials. 

SEBASS  performed  well  during  the  exercise.  Instrument  characterization  and 
calibration  efforts  allowed  spectral  analysis  to  be  conducted  using  ground  truth  spectra 
collected  by  other  instruments  at  different  times.  Although  several  instrumental  artifacts 
were  discovered  during  the  analysis  of  the  WESTERN  RAINBOW  data,  they  did  not 
severely  degrade  the  sensor’s  ability  to  collect  meaningful  data.  These  artifacts  were 
dealt  with  in  software  corrections.  The  proposed  solution  for  the  amplifier  gain  variation 
worked  well.  The  thermal  drift  and  nonlinear  response  of  the  focal  plane  were  effectively 
characterized  through  the  use  of  multiple  calibration  measurements. 

The  plastic  raler  atmospheric  compensation  method  was  designed  to  operate  on 
data  which  includes  a  large  amount  of  near-blackbody  materials.  The  sparse  vegetation 
within  the  WESTERN  RAINBOW  collection  sites  significantly  impacted  the  algorithm’s 
results.  The  small  number  of  near-blackbody  pixels  caused  the  projection  of  spectral 
features  from  the  majority  elements  in  the  scene,  such  as  sand,  into  the  derived 
atmospheric  parameters.  The  parameters  then  passed  these  features  onto  the  apparent 
emissivity  spectra.  In  addition,  the  author  found  that  the  routine  was  very  sensitive  to 
changes  in  the  number  of  bins  used  for  the  statistical  determination  of  the  transmittance 
and  upwelling  radiance  parameters.  Changes  in  binsize  as  small  as  5  or  10  pixels 
produced  changes  in  the  fine  structure  of  the  resulting  atmospheric  transmittance  and 
upwelling  radiance  profiles.  Again,  these  errors  were  propagated  through  to  the  derived 
emissivity  spectra  as  the  data  were  converted  to  radiance  at  the  ground. 

Principal  component  analysis  demonstrated  that  SEBASS  was  capable  of 
collecting  thermal  spectral  data  with  sufficient  spectral  and  radiometric  resolution  to 
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detect  and  discriminate  military  targets.  The  analysis  revealed  that  the  data  collected  by 
SEB  ASS  were  not  only  capable  of  discriminating  between  the  majority  background  and 
man-made  materials,  but  were  also  capable  of  discriminating  between  classes  of  man¬ 
made  materials.  The  analysis  showed  a  modest  increase  in  noise  introduced  by  the 
conversion  from  radiance  to  apparent  emissivity. 

The  spectral  angle  mapper  and  spectral  matched  filter  demonstrated  that  the 
conversion  to  apparent  emissivity  produced  accurate  classification  results  when  ground 
tmth  spectra  were  used  as  reference  spectra.  This  was  an  important  result  because 
atmospheric  compensation  and  the  derivation  of  apparent  emissivity  introduced 
additional  artifacts  which  may  have  obscured  or  altered  characteristic  spectral  features. 
The  classification  algorithms  demonstrated  the  successful  detection  of  various 
camouflage  materials.  They  also  were  successful  in  categorizing  the  major  terrain 
features  of  the  desert  including  vegetation,  desert  varnish,  disturbed  sand,  and  roads  and 
trails.  Of  the  two,  the  spectral  matched  filter  produced  the  best  results.  The  spectral 
angle  mapper  had  difficulty  distinguishing  between  several  of  the  camouflage  materials 
and  the  T-72  target  spectra.  It  also  had  difficulty  correctly  classifying  the  lighter  desert 
sand  and  vegetation.  It  is  believed  that  overestimation  of  surface  temperatures  resulted  in 
slope  variations  in  the  derived  emissivity  spectra.  In  the  case  of  the  above  materials, 
these  variations  produced  variations  in  spectral  reference  angle  and  resulted  in 
misclassification. 

The  spectral  matched  filter  performed  better  than  the  spectral  angle  mapper.  The 
process  used  to  convert  radiance  units  to  apparent  emissivity  often  overestimated  surface 
temperature.  This  caused  slope  variations  in  derived  emissivity  spectra.  It  is  believed 
that  the  spectral  angle  mapper  misclassifications  were  caused  by  these  slope  variations 
which  were  larger  than  the  characteristic  spectral  features  of  many  of  the  materials  within 
the  a  scene. 

The  capability  demonstrated  by  SEB  ASS  during  Exercise  WESTERN  RAINBOW 
has  many  implications  to  military  remote  sensing  at  a  time  when  battlefield  awareness  has 
become  a  dominant  factor  in  warfare.  The  results  of  this  analysis  indicate  that  an  infrared 
spectral  sensor  can  satisfy  several  key  military  intelligence  needs  such  as  target  detection 
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and  terrain  categorization.  A  thermal  spectral  sensor  like  SEB  ASS  can  detect  and 
discriminate  between  targets  under  all  lighting  conditions  based  on  spectral  differences 
alone.  It  can  also  serve  as  a  cueing  sensor  performing  anomaly  detection  over  a  large 
area.  Finally,  infrared  spectral  information  has  the  capability  to  identify  specific  target 
materials. 
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APPENDIX  A.  FIGURES 
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Figure  2.1.  US  and  Soviet,  Camouflage  and  Healthy  Green  Vegetation  Comparisons. 

From  Multispectral  Users  Guide,  1995. 
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Figure  2.3.  Classes  of  Spectral  Imaging  Sensors.  From  Multispectral  Users  Guide,  1995. 
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Figure  2.4.  Planck  Radiation  Function. 


Figure  2.5.  Reflectance  of  Pure  Water.  From  Anderson,  et  al,  1994 


Figure  2.6.  Emissivity  of  Several  Minerals.  From  Zissis, 
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Figure  2.7.  Reflectance  Spectra  of  Healthy  and  Senescent  Grass.  From  John  Hopkins 

University,  1995. 


Figure  2.8.  Sources  of  LWIR  Radiation.  After  Multispectral  Users  Guide,  1995. 


Short  wavelength  infrared 


Figure  2.9.  Atmospheric  Transmission  Characteristics.  From  Elachi,  1987. 


Figure  2.10.  Empirical  Line  Method.  After  Farrand,  1994. 


82 


Band  1 

Figure  2.14.  Illustration  of  Spectral  Angle.  After  Kruse,  et  al,  1993. 


BRNO  2 


Figure  2. 15.  Orientation  of  Signal  Vector,  b  ,  and  Clutter  Density  in  Two-Dimensional 

Band  Space.  From  Stocker,  et  al,  1990. 
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COMPONENT  1 


Figure  2.16.  Orientation  of  Signal  Vector,  b  ,  and  Clutter  Density  Principal  Component 

Space.  From  Stocker,  et  al,  1990. 
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Figure  2.17.  Weighted  Combination  of  Components  from  Figure  2.16  Which  Maximizes 
Signal-to-Clutter  Ratio.  From  Stocker,  et  al,  1990. 


Figure  3.1.  TIMS  Image  of  Cuprite  Mining  District  Processed  with  Decorrelation 
Stretch.  Abbreviations:  B  basalt,  C  carbonates,  O  opalized  units,  Q  small  blue  fan,  S 
silicified  rocks,  SS  harkless  siltstone,  T  unaltered  tuffs.  From  Kahle  and  Goetz,  1983 
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Figure  3.2.  Optical  Layout  of  BASS.  From  Bongiovi,  Hackwell,  and  Hayhurst,  1995. 


Figure  3.3.  BASS  Twin  Otter  Configuration.  From  Bongiovi,  Hackwell,  and  Hayhurst, 

1995. 
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Figure  3.5.  Spectral  Templates  of  Target  Materials.  From  Bongiovi,  Hackwell,  and 

Hayhurst,  1995. 


Figure  3.6.  Response  Plot  for  Initial  Training  Scene  of  Military  Vehicles.  From 
Bongiovi,  Hackwell,  and  Hayhurst,  1995. 
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Figure  3.7.  Response  Plot  for  Tank  Concealed  in  Trees.  From  Bongiovi,  Hackwell  and 

Hayhurst,  1995. 


Figure  3.8.  Optical  Layout  of  SEB ASS.  From  Hackwell,  1995. 
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Figure  3.9.  Diagram  of  SEBASS  Pushbroom  Imaging.  From  Hackwell,  1995. 
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Figure  3.11.  Transmittance  Spectrum  of  Polymer  Film  Used  for  Calibration.  From 

Bongiovi,  1996. 
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Figure  4.1.  Diagram  of  Malpai  Collection  Site.  From  Boisjolie,  1996. 
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Figure  4.2.  Diagram  of  JCQD  Collection  Site.  From  Boisjolie,  1996. 
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Figure  5.1.  Photo  of  JCCD  Collection  Site.  From  WESTERN  RAINBOW  Photo  CD," 

1996. 
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Figure  5.6.  SEBASS  Data  Preparation  Process. 
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Figure  5.8.  Results  of  Plastic  Ruler  Algorithm  (Without  Screening)  on  SEE  ASS  Data 

Collected  Over  Malpai  Site. 
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Figure  5.1 1.  Comparison  of  JHU  Ground  Truth  and  SEBASS-derived  Emissive  Spectra. 
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Figure  6. 1 .  Simulated  FLIR  image  of  Malpai  Collection  Site,  1000  Local  Time, 

20  October  1995. 
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Figure  6.2.  Surface  Temperature  Image  of  Malpai 

20  October  1995 
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Figure  6.3. Eigenvalues  of  Radiometrically  Calibrated  Data  and  Atmospherically 
Corrected  Data  from  Malpai  Collection  Site,  1000  Local  Time,  20  October  1995. 


Figure  6.4.  Percentage  of  Variance  Contained  in  Principal  Components  of 
Radiometrically  Calibrated  Data  and  Atmospherically  Corrected  Data  from  Malpai 
Collection  Site,  1000  Local  Time,  20  October  1995. 
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Figure  6.5.  Comparison  of  Simulated  FLIR  Image  with  First  Principal  Component  Image 
of  Radlometrically  Calibrated  Data  and  First  Principal  Component  Image  of 
Atmospherically  Corrected  Data  from  Malpai  Collection  Site,  1000  Local  Time,  20 

October  1995. 
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Figure  6.6.  Comparison  Between  Fifth  Principal  Component  Image  of  Radiometrically 
Calibrated  Data  and  Fifth  Principal  Component  Image  of  Atmospherically  Corrected  Data 
from  Malpai  Collection  Site,  1000  Local  Time,  20  October  1995. 
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Figure  6.7.  Comparison  Between  Simulated  FLIR  Image  and  Second  Principal 
Component  Image  of  Radiometrically  Calibrated  Data  from  Malpai  Collection  Site,  1000 

Local  Time,  20  October  1995. 
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Figure  6.8.  Third  Principal  Component  Image  of  Radiometrically  Calibrated  Data  from 
Malpai  Collection  Site,  1000  Local  Time,  20  October  1995. 
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Figure  6.1 1.  Red-Green-Blue  Principal  Component  Image  of  Radiometrically  Calibrated 
Data  from  Malpai  Collection  Site,  1000  Local  Time,  20  October  1995. 
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Figure  6.13.  Weighting  of  Spectral  Bands  to  Generate  Principal  Components  One,  Two, 
and  Five  of  Malpai  Collection  Site,  1000  Local  Time,  20  October  1995. 
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Figure  6.14.  First  Principal  Component  Image  of  Emissivity  Data  from  Malpai  Collection 

Site,  1000  Local  Time,  20  October  1995. 
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Figure  6.15.  Fifth  Principal  Component  Image  of  Emissivity  Data  from  Malpai  Collection 

Site,  1000  Local  Time,  20  October  1995. 
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Figure  6.16.  Red-Green-Blue  Principal  Component  Image  of  Emissivity  Data  from 
Malpai  Collection  Site,  1000  Local  Time,  20  October  1995. 
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Figure  6.17.  Results  of  SAM  Classification  of  Malpai  Collection  Site,  1000  Local  Time, 
20  October  1995  Using  Ground  Truth  Spectra. 


Figure  6.18.  Results  of  SAM  Classification  of  Malpai  Collection  Site,  1000  Local  Time, 
20  October  1995  Using  Ground  Truth  and  Extracted  Spectra. 
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Figure  6.19.  Results  of  SMF  on  Malpai  Data,  1000  Local  Time,  20  October  1995  Using 
T-72  and  T-80  Spectra.  Dark  Areas  Represent  Matches  to  Reference  Spectrum. 
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Figure  6.20.  Results  of  SMF  on  Malpai  Data,  1000  Local  Time,  20  October  1995  Using 
,  Desert  Varnish  Spectrum.  Bright  Areas  Represent  Matches  to  Reference  Spectrum 
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Figure  6.22.  Simulated  FLIR  image  of  Decoy  Line  One  of  JCCD  Collection  Site,  1000 

Local  Time,  23  October  1995. 
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Figure  6.23.  Comparison  of  FLIR  Image  of  Decoy  One  with  FT  JR  image  of  T-72  from 

the  Malpai  Data. 


Figure  6.24.  Second  Principal  Component  of  Radiometrically  Calibrated  Data  from 
JCCD  Decoy  Line  One,  1000  Local  Time,  23  October  1995. 
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Figure  6.25.  Red-Green-Blue  Principal  Component  Image  of  Radiometrically  Calibrated 
Data  from  JCCD  Decoy  Line  One,  1000  Local  Time,  23  October  1995. 
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Figure  6.26.  Results  of  SAM  Classification  of  JCCD  Decoy  Line  One,  1000  Local  Time, 
20  October  1995  Using  Ground  Truth  and  Extracted  Spectra. 
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Figure  6.27.  Results  of  SMF  on  JCCD  Decoy  Line  One,  1000  Local  Time,  20  October 
1995  Using  LOSS  Woodland  Camouflage.  Dark  Areas  Represent  Matches  to  Reference 

Spectrum. 
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Figure  6.29.  Second  Principal  Component  of  Radiometrically  Calibrated  Data  from 
JCCD  Decoy  Line  Two,  1000  Local  Time,  23  October  1995. 
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Figure  6.30.  Red-Green-Blue  Principal  Component  Image  of  Radiometrically  Calibrated 
Data  from  JCCD  Decoy  Line  Two,  1000  Local  Time,  23  October  1995. 
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Figure  6.31.  Red-Green-Blue  Principal  Component  Image  of  Emissivity  Data  from  JCCD 
Decoy  Line  Two,  1000  Local  Time,  23  October  1995. 
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Figure  6.32.  Results  of  SAM  Classification  of  JCCD  Decoy  Line  Two,  1000  Local  Time, 
23  October  1995  Using  Ground  Truth  and  Extracted  Spectra. 
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Figure  6.33.  Results  of  SAM  Classification  of  JCCD  Decoy  Line  Two  Following 
Application  of  Spatial  Filtering  Using  Ground  Truth  and  Extracted  Spectra. 
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Figure  6.36.  Red-Green-Blue  Principal  Component  Image  of  Radiometrically  Calibrated 
Data  from  JCCD  Vehicle  Target  Line,  1000  Local  Time,  23  October  1995. 
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Figure  6.37.  Results  of  SAM  Classification  of  JCCD  Vehicle  Target  Line,  1000  Local 
Time,  23  October  1995  Using  Ground  Truth  and  Extracted  Spectra. 


Figure  6.38.  Results  of  SMF  on  JCCD  Vehicle  Target  Line,  1000  Local  Time,  20 
October  1995  Using  Air  Force  Desert  Camouflage.  Dark  Areas  Represent  Matches  to 

Reference  Spectrum. 
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Figure  6.39.  Results  of  ARTSPEC  Fitting  Routine  for  SEE  ASS  Data  Collected  Over 
Arroyo  Collection  Site,  1000  Local  Time,  23  October  1995. 
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Figure  6.41.  Second  and  Third  Principal  Components  of  Radiometrically  Calibrated  Data 
from  Arroyo  Collection  Site,  1000  Local  Time,  23  October  1995. 
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Figure  6.42.  Results  of  SAM  Classification  of  Arroyo  Collection  Site,  1000  Local  Time, 
23  October  1995  Using  Ground  Truth  and  Extracted  Spectra. 
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Reference  Spectrum  Reference  Spectrum 


Figure  6.43.  Results  of  SMF  on  Arroyo  Collection  Site,  1000  Local  Time,  20  October 
1995  Using  LCSS  Woodland  Camouflage  and  LCSS  Desert  Camouflage.  Dark  Areas 
Represent  Matches  to  Reference  Spectrum. 
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APPENDIX  B.  ATMOSPHERIC  COMPENSATION  ALGORITHM 


;  PROGRAM;  plastic_ruler.pro 
;  Revision:  21  May  96 

;  PERFORMS  PLASTIC  RULER  ATMOSPHERIC  COMPENSATION 

;  Performs  basic  plastic  ruler  routine  then  repeats  using 

;  lambda  at  maximum  of  transmittance  for  reference  wavelength. 

;  Based  on  IDL  program  written  by  Bongiovi,  Hackwell,  and  Hayhurst 

;  Modified  on  21  May  to  skip  reference  calculation. 

;  FUNCTIONS  CALLED: 

;  wave.pro  FUNCTION  TO  CALCULATE  WAVELENGTH  FOR  EACH  PIXEL 

'■>  SEBASS  ARRAY  (from  Aerospace  Corporation) 

;  blackbody.pro  FUNCTION  TO  CALCULATE  BLACKBODY  RADIANCE  FOR 

;  EACH  PDCEL  IN  ARRAY 

;  brightness.pro  FUNCTION  TO  GENERATE  BRIGHTNESS  TEMP  FOR 

;  RADIANCE(uw/sr  um  cm''2) 

;  noise.estimate.pro ;  FUNCTION  TO  PERFORM  NOISE  ESTIMATION 

;  REQUIRED  FOR  PLASTIC  RULER 

;  ATMOSPHERIC  COMPENSATION 

;  screen.pro  FUNCTION  TO  SCREEN  DATA  FOR  SPECIFIC  THRESHOLD 

;  OF  BRIGHTNESS  VALUES 

;  do_plot.pro  PROCEDURE  TO  PLOT  STRIPS  OF  ATMOSPHERIC  WINDOW 

;  ksjine.pro  PROCEDURE  WRITTEN  BY  JOHN  HACKWELL  TO  PERFORM 

;  LINE  FIT  FOR  CLUSTER  OF  POINTS 

INPUT: 

;  #####_cal.dat  FLTARR( 128, 128,2000)  Calibrated  data  file 

;  cal####.cal.bip  FLTARR(128,128,256)  Calibrated  calibration  file 

;  OUTPUT: 

;  atmosphere.dat  [transmittance(128),upwelling(128)]  FLOAT 

;  results  of  basic  Hackwell 

;  DEFINE  WINDOW  AND  SET-UP  FILES 

calfile=Vrsl_02/collins/b_tmp/23oct95_1000/cal2_15.75.car 

infile=7rsl_02/collins^_tmp/23oct95_1000/jccd4/jccd4.cal’ 

outfile=Vtmp/jccd4_atmosphere.dat’ 

sand_threshold=.3 

vegetation_threshold=.5 

print,’Running  Plastic  Ruler  with: 

print,’  sand_threshold=  ’,sand_threshold 

print,’  vegetation_threshold=  ’,vegetation_threshold 

spatial=[0,127] 

frame=[200,1200] 
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spawin=indgen(spatial(  1  )-spatial(0)+ 1  )+spatial(0) 
framewin=indgen(frame(  1  )-frame(0)+ 1  )+fTame(0) 

print,  ’Window  Contains  n_elements(spawin)*n_elements(framewin),’  pixels’ 

;  SETUP  WAVELENGTH  ARRAY 

wv=wave() 

lambda=fltarr(  1 28,n_elenients(spawin)) 

for  i=0,n_elements(spawin)-l  do  lambda(*,i)=wv(*,63) 

;  PERFORM  NOISE  ESTIMATE 

noise=noise_estimate(calfile, spatial) 

;noise2=noise_estimate(calfile2, spatial) 
plot, noise 

;oplot,noise2,linestyle=2 

stop 

;  LOAD  IN  WINDOW  OF  DATA 

openr,  unit,  infile, /getjun 

temp=fltarr(  1 28, 1 28) 
for  loop=0,  frame(0)-l  do  begin 
forrd,unit,temp 

endfor 

data=fltarr(  1 28, 1 28,n_elements(framewin)) 

forrd,unit,data 
close,  unit 
ffeejun,unit 

atmdata=  data(*,spatial(0):spatial(l),*) 
allatmdata=reform(atmdata,128,n_elements(spawin)*$ 
n_elements(framewin)) 

;  SCREEN  FOR  BLACKBODIES 

print, ’screening  for  blackbodies’ 

;  CHECK  FOR  SAND  MIXING  (10  to  9.45) 

;  average  pixels  ivo  10  and  compare  to  pixels  ivo  of  9.45 

;  which  is  in  reststrahlen  band 

lambdaA=[9.9,10.1] 

lambdaB=[9.3,9.6] 

vegetation=screen(allatmdata,lambda,lambdaA,lambdaB,sand„threshold) 

vegetationcount=n_elements(vegetation) 

print,’  Selected  ’,vegetationcount,’  vegetation  pixels’ 

;  CHECK  FOR  SENESCENT  GRASS  (10.65  to  1 1.5) 

;  average  pixels  ivo  10.65  and  compare  to  pixels  ivo  of  1 1.5 

;  which  is  in  peak  in  senescent  grass  (from  JHU) 
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lambdaA=[l  0,6,10.7] 
lambdaB=[l  1.4,1 1.6] 

blackbody_veg=screen(allatmdata(*, vegetation), lambda, lambdaA,lambdaB,  $ 
vegetation_threshold) 
selected=vegetation(blackbody_veg) 
count=n_elements(selected) 
print,’  Selected  count,  ’  pixels’ 

;  DISPLAY  PLOT  OF  STRIP 

one_two_three=’  1  ’ 

read,  ’Select  a  plot  option:  (1)  none,  (2)  postscript,  (3)  screen:  ’,$ 
one_two_three 


case  1  of 

(one_two_three  eq  2):  begin 
set_plot,’ps’ 

device, file=’window.ps’,/helvetica,/boid,bits=8 
do_plot,infile,refomi(atmdata(63,*,*)),spawin,framewin,  $ 
2,25,vegetation,’ All  Vegetation’ , selected, ’Healthy ’ 
device,/close_file 
set_plot,’x’ 
end 

(one_two_three  eq  3):  begin 
set__plot,’x’ 

device, true=24,retain=2 
zoom=950/  n_elements(framewin) 
if  (zoom  eq  0)  then  zoom=l 
window, 0,xsize=3*zoom*n_elements(spawin)+20,  $ 
ysize=zoom*n_elements(ffamewin)+50 
do_plot,infile,reform(atmdata(63,*,*)),spawin,framewin,  $ 
zoom,  1 ,  vegetation,  ’Vegetation’ ,  selected,  ’  Screened  ’ 
end 

else:  print,” 
endcase 

read, ’Enter  Binsize:  ’,binsize 

;  FINDTMAX 

bright=fltarr(  1 28, count) 
tbright=fltarr(count) 

locations=intarr(count) 

for  i=0, count- 1  do  begin 

bright(*,i)=brightness(alIatmdata(*,selected(i)),lambda(*,0)) 

tbright(i)=max(rebin(bright(*,i),64)) 

shortlambda=rebin(lambda,64) 

locations(i)=shortlambda(!c) 

endfor 

histo=histogram(locations,min=:0) 

histomax=max(histo) 
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print,’  Most  Tmax  at  ’,shortlambda(!c,0) 

;  BASIC  PLASTIC  RULER  FIT! 
print,  ’Doing  BASIC  correction’ 

;  MAKE  BLACK  BODIES 
allblackbodies=fltarr(  1 28,  count) 
for  i=0,count- 1  do  begin 

allblackbodies(*,i)  =  blackbody(tbright(i),lambda(*,0)) 

endfor 


DO  CORRECTION  CALC 

HACKWELL:  temp=fltarr(5,128) 
correction=fltarr(5, 128, 128) 

for  i=0,127  do  begin 
temp(*,i)=  $ 

ks__line(allblackbodies(i,indgen(count)),  $ 
allatmdata(i, selected), noise(i),binsize,f) 

endfor 

transmittance=temp(  1  ,*) 
upwelling=temp(0,*) 

maxtrans=max(transmittance) 

channel=!c 

print,’  Max  transmittance  at  ’,  wv(!c,63) 

;  PLOT  TRANSMITTANCE  CURVES 

window,! 

!p.multi=:[0  ,1,2] 

plot, wv(*,63), transmittance, title=’ Atmospheric  Transmittance’  ,$ 
ytitle=’  Transmittance’  ,xtitle=’ Wavelength  (microns)  ’ 
xyouts,  8,1.1, infile 

text=’Screening  Thresholds:  ’  +  string(sand_threshold)  +  $ 
string(vegetation_threshold) 
xyouts,8,.25,text 
text=’Binsize:’  +  string(binsize) 
xyouts, 8,.15, text 

plot, wv(*,63),upwelling,title=’ Atmospheric  Upwelling  Radiance’, $ 
ytitle=’Upwelling  Radiance’, xtitle=’ Wavelength  (microns)’ 

window, 2 
!p.multi=0 

plot,wv(*,63),(allatmdata(*,selected(100))-upwelling)/transmittance,/ynoz 

oplot,wv(*,63),(allatmdata(*,selected(200))-upwelling)/transmittance 

oplot,wv(*,63),(allatmdata(*,selected(300))-upwelling)/transmittance 

;  PERFORM  LOOP 
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read, ’Select  new  binsize  or  zero  to  stop:  input 
if  (input  ne  0)  then  begin 
binsize=input 
goto,  HACKWELL 

endif 

;  SAVE  RESULTS 
one_or_two=’ r 


openw,unit,  outfile,/get_lun 

print, ’Saving  BASIC  results’ 

forwrt,unit,transmittance 

forwrt,unit,upwelling 

close, unit 
free_lun,unit 

end 

;FUNCTIONS  AND  PROCEDURES  CALLED  BY  ATMOSPHERIC  COMPENSATION  ALGORITHM 
;  FUNCTION:  wave.pro 

;  Generated  by  R.  Bongiovi,  Aerospace  Corporation 

;  Returns  an  array  containing  the 

;  wavelength  for  the  SEBASS  LWIR  array  using  the 

;  plastic  filter  data  from  20  Oct  95  (file:  call_48.36p.dat) 

function  wave 

;  Set  up  and  return  function 

wavel  =  fltarr(128,128) 
x  =  findgen(128) 

aO  =  57.0758  -  0.0216030*x  +  0.000100471  *x^2 
al  =  1.03258  +  1.88537e-5*x  +  9.2521  le-8*x^2 
a2  =  replicate(-0.000100671,128) 
for  j  =  0,127  do  begin 

wavel(*,j)  =  sqrt(a0(j)  +  al(j)*x  +  a2(j)*x^2) 
endfor 

return, wavel 
end 


;  FUNCTION:  blackbody.pro 
;  revision:  30  Mar  96 

;  FUNCTION  TO  GENERATE  BLACKBODY  RADIANCE  FOR  AN  INPUT  TEMP  (KELVIN) 

;  INPUT: 

;  temperature  (kelvin) 
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float  data 


wv  (um) 

float  array  stored  as  dim(l),dim(2) 
contains  wavelength  bins  for  calculation 


;  OUTPUT: 

;  radiance  (uw/um  cm^2  sr) 

;  double  float  data 

function  blackbody,temp,wv 

dim=size(wv) 

radiance=dblaiT(dim(  1  ),dini(2)) 

cl=L19104el0 

c2=L438769e4 

radiance=c  l/((wv^5)*(exp(c2/(w  v*temp))- 1 )) 

return, radiance 
end 

;  FUNCTION:  brightness.pro 
;  revision:  15  AUg  96 

;  FUNCTION  TO  GENERATE  BRIGHTNESS  TEMP  FOR  RADIANCE(uw/sr  um  cm^2) 

;  INPUT: 

;  radiance  FLOAT  (,)  Radiance  at  each  pixel 

;  wav  FLTARR  (,)  Wavelengths  for  each  pixel 

;  OUTPUT: 

;  bright  FLTARR(128,128)  [kelvin] 

function  brightness, radiance, wav 
dim=size(wav) 

brightness=dblarr(dim(  1  ),dim(2)) 

cl=1.19104el0 

c2=1.438769e4 

brightness=c2/(wav*alog(c  l/(way^5  *radiance)+l  .0)) 

return, brightness 
end 

;  PROCEDURE:  screen.pro 
;  Revision:  7  May  96 

;  FUNCTION  TO  SCREEN  DATA  FOR  SPECIFIC  THRESHOLD  OF 
;  BRIGHTNESS  VALUES 

;  INPUT: 

;  data  FLTARR(  128, range)  (array  of  radiance  data) 

index  INT  (index  of  data  elements  to  check) 

;  lambda  FLTARR(SEBASS  Wavelengths) 
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lambdaA  FLTARR(2)  (range  A) 
lambdaB  FLTARR(2)  (range  B) 
threshold  FLT  (threshold  value) 


;  OUTPUT: 

;  noise  [noise(128,n_eiements(SPAWIN))]  FLTARR 

function  screen,data, lambda, lambdaA, lambdaB , threshold 

x=size(data) 

number=x(2) 

vv^v=wave() 

tempA=fltarr(number) 

tempB=fltarr(number) 

for  i=0L,  number- 1  do  begin 

x=where(wv(*,63)  gt  lambdaA(O)  and  wv(*,63)  It  lambdaA(l),xcount) 
tempA(i)=total(brightness(data(x,i),lambda(x,0)))/xcount 
x=where(wv(*,63)  gt  lambdaB(O)  and  wv(*,63)  It  iambdaB(l),xcount) 
tempB(i)=total(brightness(data(x,i),lambda(x,0)))/xcount 

endfor 

selection=>vhere(abs(tempA-tempB)  It  threshold) 

return,  selection 

end 

;  PROCEDURE:  do_plot.pro 
;  Revision:  9  May  96 

PROCEDURE  TO  PLOT  STRIPS  OF  ATMOSPHERIC  WINDOW 
;  INPUT: 

;  infile  STR  (name  of  file) 

;  strip  FLTARR(spectral,  spatial) 

;  spawin,framewin  FLTARR(spatal  and  frame  window  index) 

;  zoom  FLT  (scale  factor) 

;  pos^factor  FLT  (multiplicative  factor  for  PS  plot) 

;  screenA,screenB  FLTARR(selected  pixels) 

;  nameA,nameB  STR  (name  of  selections) 

;  OUTPUT: 

pro  do_plot,infile,strip,spawin,framewin,zoom,pos_factor,screenA,nameA,  $ 
screenB,nameB 

pIotstrip=rebin(strip,n_elements(spawin)*zoom,n_elements(framewin)*zoom) 

min=min(plotstrip) 

max=max(piotstrip)+.2*median(plotstrip) 

plotstrip(0,0)=max 

tvscl,bytscl(plotstrip,min,max),0,40*pos_factor 

alLplotstrip=reform(strip,n_elements(spawin)*n_elements(framewin)) 

all_plotstrip(screenA)=max 
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plotstrip=refonn(alLplotstrip,n_elements(spawin),n_elements(framewin)) 
plotstrip=rebin(plotstrip,n_elements(spawin)*zoom,n_eiements(framewin)*zoom) 
tvscl,bytscl(plotstrip,min,max),zoom*pos_factor*n_elements(spawin)+10,  $ 

40*pos_factor 

alLplotstrip=reform(strip,n_elements(spawin)*n_elements(£ramewin)) 

alLplotstrip(screenB)=:max 

plotstrip=reform(all_plotstrip,n_elements(spawin),n_elements(fTamewin)) 
plotstrip=rebin(plotstrip,n_elements(spawin)*zoom,n_elements(framewin)*zoom) 
tvscl,bytscl(plotstrip,min,max),2*zoom*pos_factor*n_elements(spawin)+20,  $ 

40*pos_factor 

xyouts, 0.15, .04, ’Window’, alignment  =  0.5, size  =  L5,/normal 
xyouts,0.5,. 04, nameA, alignment  =  0.5, size  =  1.5, /normal 
xyouts,0.85, . 04 ,nameB, alignment  =  0.5, size  =  1.5, /normal 
xyouts,.5,.01,infile,alignment  =  0.5,size=.75  ,/normal 

end 

;  FUNCTION:  noise_estimate.pro 
;  Revision:  31  Mar  96 

;  FUNCTION  TO  PERFORM  NOISE  ESTIMATION  REQUIRED  FOR  PLASTIC  RULER 
;  ATMOSPHERIC  COMPENSATION 

;  INPUT: 

;  filename  STR(name  of  calibration  file) 

;  spatial  intaiT(2-  start  and  stop  spatial  windows) 

;  OUTPUT: 

;  noise  [noise(128,n_elements(SPAWIN))]  FLTARR 

function  noise_estimate,filename,spatial 
;  LOAD  IN  CALIBRATION  FILE 
nrecs=128 

cal_data=fltarr(  128,1 28,nrecs) 

openr,  unit,  filename,/get_lun 

forrd,unit,cal_data 

close, unit 

;  PERFORM  NOISE  ESTIMATE 

noise=fltarr(128) 

for  i=0,127  do  begin 

noise(i)=stdev(caLdata(i,spatial(0):spatial(l),*)) 

endfor 

retum,noise 

end 
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